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Abstract 

We reduce the problem of finding an augmenting path in a general graph to 
a reachability problem in a directed bipartite graph. A slight modification 
of depth-hrst search leads to an algorithm for finding such paths. Although 
this setting is equivalent to the traditional terminology of blossoms due to 
Edmonds, there are some advantages. Mainly, this point of view enables 
the description of algorithms for the solution of matching problems without 
explicit analysis of blossoms, nested blossoms, and so on. Exemplary, we 
describe an efficient realization of the Hopcroft-Karp approach for the com¬ 
putation of a maximum cardinality matching in general graphs and a vari¬ 
ant of Edmonds’ primal-dual algorithm for the maximum weighted matching 
problem. 


1 Introduction and motivation 

Since Berge’s theorem in 1957 [1] it has been well known that for constructing a 
maximum matching, it suffices to search for augmenting paths. But until 1965, 
only exponential algorithms for finding a maximum cardinality matching in non- 
bipartite graphs have been known. The reason was that one did not know how to 
treat odd cycles, the so-called “blossoms”, in alternating paths. In his pioneering 
work, Edmonds na solved this problem by shrinking these odd cycles. Because 
each shrinking reduces the number of nodes in the current graph at least by two, the 
total number of shrinkage’s is bounded by | where n is the number of nodes in the 
graphs. Hence, Edmonds original algorithm uses only O(n^) time. In [21 [T7ll?nIS7] . 
it is shown how to avoid explicit shrinking of odd cycles. All these algorithms use 
O(n^) or 0{nma{m + n,n)) time where m is the number of edges in the graph and 
a is the functional inverse of Ackermann’s function. 
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In 1973, Hopcroft and Karp [2l] proved the following fact. If one computes in 
one phase a maximal set of pairwise disjoint shortest augmenting paths and aug¬ 
ments these paths then 0{^/n) such phases would be sufficient. For the bipartite 
case they showed that a phase can be implemented using a breadth-first search fol¬ 
lowed by a depth-first search. This led to an 0{n + m) implementation of one phase 
and hence, to an 0{y/nm) algorithm for the computation of a maximum matching 
in bipartite graphs. The implementation of a phase for non-bipartite graphs is 
much harder. In 1975, Even and Kariv [IIII2S] presented an 0(min{n^, mlogn}) 
implementation of a phase leading to an 0(min{n^'®, ydim logn}) algorithm for 
the computation of a maximum matching in general graphs. Galil called the 
full paper [25] “a strong contender for the ACM Longest Paper Award”. Tarjan 
[3^ called their paper “a remarkable tour-de-force”. In 1978, Bartnik |3] gave an 
alternative O(n^) implementation in his unpublished Ph.D. thesis (see [23|)- In 
1980, Micali and Vijay Vazirani |3T] have presented an 0{m -|- n) implementation 
of a phase without the presentation of a correctness proof. The algorithm as pre¬ 
sented in [31] is not correct since their definition of “tenacity” does not work in 
all cases. In 1988, Peterson and Loui have given an informal correctness proof of 
the incorrect algorithm of Micali and Vazirani. Fourteen years later, Vijay Vazi¬ 
rani [3H] has repaired the mistake by changing the definition of “tenacity” and 
provided a correctness proof. Recently [36], Vazirani has presented a new version 
of a correctness proof. With respect to his proof in [35], he writes: “Although 
the statements of these theorems were largely correct, their proofs, which involved 
low level arguments about individual paths and their complicated intersections 
with other structures, were, in retrospect, incorrect.” In 1999 [H], I have tried 
to combine ideas of Micali and Vazirani and my framework. As pointed out by 
Ross McConnell [29], the algorithm, as described in [H], is not correct. In 1991, 
Cabow and Tarjan pni have given an efficient scaling algorithm for general graph¬ 
matching problems. As pointed out in their paper, by a slight modification of their 
algorithm, they have obtained another 0{m + n) implementation of a phase. The 
history of efficient implementations of a phase of the Hopcroft-Karp approach for 
general graphs illustrates the need of a framework which allows a clear description 
and an elaborated correctness proof of matching algorithms. 

The hrst polynomial algorithm for the maximum weighted matching problem 
has been also given by Edmonds [13]. A straightforward implementation of his 
algorithm has run time 0{n‘^m). Cabow [TB] and Lawler [27] have developed 0{n^) 
implementations of Edmonds’ algorithm. Calil, Micali and Cabow [22] have given 
an 0{nm\ogn) implementation. At the first SODA, Cabow [18] has presented an 
implementation of Edmonds algorithm which uses complicated data structures and 
stated that the time complexity of his implementation is 0{n{m -|- nlogn)). 

Our goal is to avoid sophisticated explicit analysis of (nested) blossoms. Hence, 
we reduce the problem of Ending an augmenting path to a reachability problem 
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in a directed, bipartite graph. We show, how to solve this reachability problem 
by a modihed depth-hrst search. This approach yields an algorithm which is not 
fundamentally different from previous algorithms which use Edmonds’ traditional 
terminology of blossoms. But if Edmonds’ algorithm is used as a subroutine with 
respect to the solution of a more involved problem, using the framework of the 
reachability problem which avoids the explicit consideration of blossoms can sim¬ 
plify the situation considerably. To illustrate this, we describe a realization of 
the Hopcroft-Karp approach [21] for the computation of a maximum cardinality 
matching in general graphs. Furthermore, we show how to use the modihed depth- 
hrst search algorithm in the primal step of Edmonds’ maximum weighted matching 
algorithm. For the description of the primal-dual method, we use no linear pro¬ 
gram and no duality theory for linear programs. A straightforward 0{nm\ogn) 
implementation will be described as well. 

In Section 2, the basic algorithm is presented. After the description of the 
reduction to a reachability problem in a directed, bipartite graph, this reachabil¬ 
ity problem is solved by a modihcation of depth-hrst search. The correctness of 
the algorithm is proved and an efficient implementation is given. In Section 3, 
a realization of the Hopcroft-Karp approach for general graphs is described and 
its correctness is proved using the framework of the reachability problem. Fur¬ 
thermore, an efficient implementation is given. In Section 4, we show how to use 
the modihed depth-hrst search as a subroutine in Edmonds’ maximum weighted 
matching algorithm. Furthermore, we describe an 0{nm\ogn) implementation of 
this approach. 


2 The basic algorithm 

After the reduction of the problem of hnding an augmenting path to a reachability 
problem, we shall describe the solution of the reachability problem by a modih¬ 
cation of depth-hrst search. The resulting algorithm is equivalent to Edmonds’ 
algorithm since to each run of Edmonds’ algorithm there corresponds a run of the 
modihed depth-hrst search and vice versa [33|. Hence, the correctness of the modi¬ 
hed depth-hrst search follows directly from the correctness of Edmonds’ algorithm. 
Nevertheless, we shall prove the correctness of the algorithm directly without the 
use of the correctness of Edmonds’ algorithm. Then we shall describe an efficient 
implementation of the algorithm. The description of Edmonds’ algorithm within 
the framework of the reachability problem seems not to be simpler than an elabo¬ 
rated description of Edmonds’ algorithm within the traditional framework of blos¬ 
soms. However, as we shall see later, if Edmonds’ algorithm is used as a subroutine 
with respect to the solution of a more involved problem, using the framework of 
the reachability problem can simplify the situation considerably. 
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2.1 Definitions and the general method 

A graph G = (V, E) consists of a finite, nonempty set of nodes V and a set of 
edges E. G is either directed or undirected. In the (nn-)directed case, each edge is 
an (nn-)ordered pair of distinct nodes. A graph G = (V, E) is bipartite if V can 
be partitioned into disjoint nonempty sets A and B snch that for all {u,v) G E, 
u & A and v & B or vice versa. Then we often write G = {A,B,E). A path 
P from nGldtowGldisa seqnence of nodes v = no, ni,..., wa: = w, which 
satisfies G i?, for 0 < i < A;. The length |P| of P is the nnmber k of 

edges on P. P is simple if Vi 7 ^ Vj, for 0 < i < j < k. For conveniences, P will 
denote the path no, ni,..., Vk, the set of nodes {no, ni,..., n^}, and the set of edges 
{(' 1 ^ 0 , {^ 1 , 02 ),..., (nfc-i, nfc)}. If there exists a path from n to tn (of length 1) 

then n is called a (direct) predecessor of tn, and w is called a (direct) successor of 
n. Let G = iV, E) be an nndirected graph. M C P is a matching of G if no two 
edges in M have a node in common. A matching M is maximal if there exists no 
e & E \ M snch that M U {e} is a matching. A matching M is maximum if there 
exists no matching M' C P of larger size. Given an nndirected graph G = (V, P), 
the maximum matching problem is finding a maximnm matching M O E. A path 
P = no,ni, ... ,Vk is M-alternating, if it contains alternately edges in M and in 
P \ M. A node n G G is M-free if v is not incident to any edge in M. Let 
P = ■ ■ ■ ,Vk be a simple M-alternating path. P is M-augmenting if no and 

Vk are M-free. Let P be an M-angmenting path in G. Then M ® P denotes the 
symmetric difference of M and P; i.e., M©P = M\ PUP \ M. It is easy to see 
that M © P is a matching of G, and |M © P| = |M| + 1. 

The key to most algorithms for hnding a maximnm matching in a graph is the 
following theorem of Berge |1] . 

Theorem 1 Let G = (V, P) be an undirected graph and M ^ E be a matching. 
Then M is maximum if and only if there exists no M-augmenting path in G. 

Berge’s theorem directly implies the following general method for finding a maxi¬ 
mnm matching in a graph G. 

Algorithm 1 

Input: An nndirected graph G = {V,E), and a matching M E. 

Output: A maximnm matching Mmax- 

Method: 

while there exists an M-angmenting path 

do 

constrnct snch a path P; 

M ■.= M ®P 

od; 
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The key problem is now this; How to hnd an M-augmenting path P, if such a 
path exists? We solve this key problem in the following way. 

1. We reduce the key problem to a reachability problem in a directed, bipartite 
graph Gm = {V',Em)- 

2. We solve this reachability problem constructively. 

2.2 The reduction to a reachability problem 

In the bipartite case, we construct from G = {A, B, E) and a matching M E a. 
directed graph Gm = (H U {s}, B U {t}, Em) by directing the edges in M from A 
to B, and directing the edges in i? \ M from B to A. Additionally, we add two 
new nodes s and t to AU B, add for each M-free node b G B the edge (s, b) to Em, 
and add for each M-free node a G A the edge (a, t) to Em- It is easy to prove that 
there is an M-augmenting path in G if and only if there is a simple path from s 
to t in Gm- This reachability problem can be solved by performing a depth-first 
search of Gm with start node s. Now we shall consider the general case. 

Let G = {V, E) be an undirected graph and M O E he a matching. Let 
Vm := {x G H I a; is M-free}. For the dehnition of Gm we have the following 
difficulty. Let us consider the graph described in Figure 1 where edges in M are 
wavy. The M-augmenting path from x to y enters the edge e in 6 and leaves the 
edge from a. The M-augmenting path from v to w enters e in a and leaves e from 
b. A priori, we cannot divide the set of nodes V into two sets A and B such that an 
M-augmenting path exists in G if and only if there exists an M-augmenting path, 
using alternately nodes from A and from B. Hence, for dehning Gm, we introduce 
for each node v E V two nodes [n. A] and [n, B] such that an analogous construction 
of a graph Gm is possible. Both edges ([n. A], [tc, B]) and ([tc. A], [n, B]) are in Gm 
if and only if {v,w) G M. Both edges {[x,B], [y,A\) and {\y,B], [2^)^]) are in Gm 
if and only if (x, y) E E\ M. Additionally, we add for each M-free node v E V the 
edges (s, [v,B]) and ([n,A],t) to Gm, where s and t are two new distinct nodes. 
More formally, let Gm ■= {V'^Em) where 

V' := {[v,A],[v,B] I V G H} U {s,t} s,t ^V, s 
Em ■= {{[v,A],[w,B]),{[w,A],[v,B]) \ {v,w) E M} 

hi {[x, B],[y, A]), {[y, B],[x, A]) \ {x,y) eE\M} 

U{(s, [n,H]), ([v, A],f) I V E Vm}- 

Analogously to the bipartite case, we have directed the edges in M “from A to B” 
and the edges in i? \ M “from B to A”. Since the distinct nodes [n. A] and [n, B] 
in V' correspond to the same node v in V, it does not suffice to construct a simple 
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Figure 1: Difficulty with respect to the definition of Gm- 

path from s to f in Gm for finding an M-augmenting path in G. Hence, we define 
strongly simple paths in Gm which cannot contain both nodes [n, A] and [n, B], for 
all V E V. A path P in Gm is strongly simple if 

a) P is simple, and 

b) V[n, H] e : [n. A] E P ^ [n, B] ^ P. 

Now we can formulate the reachability problem in Gm which is equivalent to the 
problem of finding an M-augmenting path in G. 

Theorem 2 Let G = (V, E) be an undirected graph, M O E be a matching, and 
Gm = {V',Em) be defined as above. Then there exists an M-augmenting path in 
G if and only if there exists a strongly simple path from s to t in Gm- 

Proof: Let Q = Wi,W 2 , - - - ,Wi_i,Wi be an M-augmenting path in G. Then 

Wi 7 ^ Wj, 1 < i < i < I, and wi,wi E Vm- Hence, by the construction of Gm, 
Q' = s, [wi,B], [w 2 , H],..., [wi-i, B], [wi, H], f is a strongly simple path in Gm- 

Let P = s, [vi,B], [^2,-4], [v 3 ,B ],..., [vk-i,B], [vk,A],t be a strongly simple 
path in Gm- Then Vi Vj, 1 < i < j < k, and vi,Vk E Vm- Hence, P' = 
Vi,V 2 , - - - ,Vk is an M-augmenting path in G. ■ 

Because of Theorem the reachability problem equivalent to the problem of 
finding an M-augmenting path in G is to find a strongly simple path from s to f 
in Gm- 
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2.3 The solution of the reachability problem 

Depth-first search (DFS) hnds simple paths in a directed graph. Hence, we cannot 
use DFS directly for the solution of the reachability problem in Gm- We will modify 
the usual DFS such that the modihed depth-hrst search (MDFS) hnds precisely the 
strongly simple paths in Gm- Let [v, H] := [v,B] and [v,B] := [v,A\. Remember 
that a DFS partitions the edges of the graph into four categories [T]. Similarly, the 
edges of Gm are partitioned into hve categories by an MDFS of Gm'- 

1. Tree edges, which are edges leading to new nodes [v,W], X G {A,B} for 
which [v,X] is not a predecessor during the search. 

2. Weak back edges, which are edges leading to new nodes [n, A\ for which [n, B] 
is a predecessor during the search. 

3. Back edges, which go from descendants to ancestors during the search. 

4. Forward edges, which go from ancestors to proper descendants but are not 
tree edges. 

5. Cross edges, which go between nodes that are neither ancestors nor descen¬ 
dants of one another during the search. 

Like DFS, MDFS uses a stack K for the organization of the search. Analogously 
to DFS, the MDFS-stack K dehnes a tree, the MDFS-tree T. Before describ¬ 
ing MDFS in detail, we shall describe the algorithm informally. TOP(iF) de¬ 
notes the last node added to the MDFS-stack K. In each step, MDFS con¬ 
siders an edge (TOP(iF), [tc, F]) which has not been considered previously. Let 
e = ([u,X], [w,^]) be the edge under consideration. We distinguish two cases. 


X = 

A, i.e. (u, w) e M. 

tree 

edge 

X = 

B, i.e. {v,w) e E\M 



2.1 

[w. A] e K 

back 

edge 

2.2 

[ic. A] ^ K but [w, B] G K 




i) [w. A] has been in K previously 

cross 

edge 


ii) [w. A] has not been in K previously 

weak back 

edge 

2.3 

[w. A] ^ K and [w, B] ^ K 




i) [w. A] has been in K previously 

forward or cross 

edge 


ii) [tc. A] has not been in K previously 

tree 

edge 
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MDFS differs from DFS only in Cases 2.2.ii and 2.3.i. Next, we shall discuss both 
of these cases. 

Case 2.2.ii: Since [w, A] has not been in K previously, DFS would perform the op¬ 
eration PUSH([t(;, A]). Since [tc, B] & K and MDFS should only construct strongly 
simple paths in Gmi the operation PUSH([t(;, A]) is not performed by MDFS. 

Note that the path P = s, Pi,[w, B], P 2 ,[v, B] dehned by the MDFS-stack 
K is strongly simple. Hence, the path Q = [x, A], P 22 ,[v, B],[w, A] where P 2 = 
P 2 i,[x, A], P 22 is strongly simple for each node [x,A\ on P 2 . We say then that 
MDFS has found the strongly simple path Q from [x. A] to [tc. A]. Since the path 
Q is above the node [w,B] with respect to the MDFS-stack K, after the execution 
of the operation POP{[w, B]), no node on Q is in the MDFS-stack K. Moreover, 
as we shall prove later, for all nodes [ziX] on P 2 the operations PUSH{\z, X]), 
POP{[z,X]), PUSH{[z,X]) and POP{[z,X]) have been performed before the op¬ 
eration POP{[w^B]). 

Case 2.3.i: Since [tc,H] has been in K before, DFS would perform no PUSH- 
operation. But the different treatment of Case 2 . 2 .ii can cause the following situa¬ 
tion: MDFS has found a strongly simple path Q = [tc, H], Q\ [m, A\ from the node 
[w,A\ to a node [u,A\ but at that moment, the node [u,B] was below [w,A\ in 
the MDFS-stack K such that the operation PUSH([m, H]) has not been performed. 
But now, [ti, B] ^ K. 

As we shall prove later, the paths P from s to [u, B] and Q from [w, A] to 
[m. A] are strongly disjoint] i.e. there is no node [r, X] on P, X E {A, H} such that 
{[r. A], [r, H]} n Q 7 ^ 0. Since MDFS has found a strongly simple path P, Q from s 
to [m. A], MDFS now performs the operation PUSH([m, A]). 

Note that with respect to depth-hrst search, the DFS-stack contains exactly 
the current search path. With respect to the modihed depth-hrst search, the sit¬ 
uation is different. In Case 2.3.i, the node [u, A] is pushed. But to obtain a cur¬ 
rent search path, between the nodes [v,B] and [m. A], we have to insert any path 
[tc, A],(5', [m. A] which has been found by MDFS. Since we do not want to forget 
the information about the hrst node [w, A] on the path which we add between the 
nodes [u, B] and [u, A], we create the artihcial tree edge ([v, B], [u, A])[^^a]. Such an 
edge is called extensible edge. It is possible that there exists various such paths Q'. 
Hence, after the execution of PUSH([m, A]), the number of corresponding current 
search paths can increase. Whenever we say that we consider a eurrent seareh path 
we mean that we can take an arbitrary corresponding current search path. If we 
add to the constructed MDFS-tree T all forward, back, cross and weak back edges 
and every extensible edge ([v, B], [u, A])[u,^^] is replaced by all strongly simple paths 
Q = [w, B], [w. A], Q', [u, A], then we obtain the expanded eurrent MDFS-tree T^xp- 

We say that MDFS has construeted a path P if the MDFS-stack K contains the 
path P where each extensible edge ([u, H], [m, A\)\^^a] in K is replaced by one of the 



strongly simple paths Q = [n, B], [w, A\,Q\ [ti, A\ which replace this extensible edge 
in Tf,xp- We say that MDFS has formed a strongly simple path P if contains P. 
We say that MDFS has found a strongly simple path P', [n, B], [tc, A] if the path 
P',[v,B] is formed by MDFS and the edge {[v, B],[w, A]) is a considered weak 
back edge. Next we shall describe MDFS in detail. We have to solve the following 
problem; How to hnd the node [u, A] in Case 2.3.i? For the solution of this problem, 
we assume that MDFS is organized such that for all nodes [w, A] G W, the following 
holds true: 

After performing the operation POP([w, A]), MDFS has always computed a set 
L[w,a] of nodes such that Li^^a] contains exactly those nodes [u, A] G V satisfying 
the requirements that 


1. MDFS has found a path P = [w, A], Q, [u, A] with [u, B] ^ Q, 

2. PUSH([m, A]) has never been performed, and 

3. POP([m,P]) has been performed. 


Before the execution of POP([tc, A]), we £x := 0- 

In the description of MDFS we assume for all [w, A] G V' that is computed 


I I ^ 


1. The 


correctly. As we shall prove later, it will always hold that 
computation of 

For [v,X] G W, N[v,X] denotes the adjacency list of [n,X]. 


Lj^A] as well as an efficient implementation of MDFS, can be 


found in Section 2.5 


Note that after the POP of the head [w,B] of a matched edge ([n. A], [w,B]), the 
POP of the tail [n. A] of this edge has also to be performed. Hence, at the end of 
the procedure SEARCH, we shall have two POP-operations. 


Algorithm 2 (MDFS) 

Input: Gm = (W,Pm)- 

Output: A strongly simple path P from s to t, if such a path exists. 

Method: 

Initialize the stack K to be empty; 

PUSH(s); 

while A 7 ^ 0 and no path from s to f is constructed 

do 

SEARCH 

od. 


SEARCH is a call of the following procedure. 


procedure SEARCH; 
if TOP(A) = t 
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then 

reconstruct a strongly simple path P from s to t 
which has been constructed by the algorithm 

else 

mark TOP(if) “pushed”; 

for all nodes [w, Y] e N[TOP{K)] 

do 

iiY = B 

(Case 1) then 

FVSR{[w,B]); 

SEARCH 

(Case 2) else 

if [tc, A] ^ K 

(Case 2.1) then 

no PUSH-operation is performed 

else 

if [w, B]eK 

(Case 2.2) then 

no PUSH-operation is performed 

(Case 2.3) else 

if [tc, A] is marked “pushed” 
(Case 2.3.i) then 

while Li^^a] 7^ 0 

do 

choose any [m, A] G 
PUSH([m,A]); 

SEARCH 

od 

(Case 2.3.ii) else 

PUSH ([w, A]); 

SEARCH 

fi 

fi 

fi 

fi 

od; 

POP; 

POP 

fi. 
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2.4 The correctness proof of MDFS 

The correctness proof of MDFS is inspired by the correctness proof of DFS. Bnt in 
contrast to DFS, the proof is difhcnlt. The difhcnlties come from the fact that the 
MDFS-stack does not contain the whole cnrrent search path and the decisions taken 
by the algorithm only depend on the content of the cnrrent stack. Hence, the proof 
that the algorithm constrncts only strongly simple paths is involved. First we shall 
prove some lemmas. The hrst lemma implies that the hrst PUSH-operation which 
destroys the property “strongly simple” mnst push a node with second component 
H. 

Lemma 1 As long as MDFS constructs only strongly simple paths, the following 
holds true: After the operation PUSH{[v, A]) where v is not M-free, the operation 
PUSH{[w, B]) where ([n, H], [w, 5]) G Em always follows without destroying the 
property “strongly simple”. 

Proof: After the execution of the operation PUSH([n, A]), always the unique edge 
([n. A], [w,B]) G Em is considered and the operation PUSH([r(;, 5]) is performed. 
If this operation would destroy the property “strongly simple”, then [w, A] and 
hence, [n, B] would be on a current search path. But then, already the operation 
PUSH([n, A]) would have destroyed the property “strongly simple”, a contradiction. 


The next lemma shows that in a certain situation MDFS constructs a path from 
s to a node [x, A]. 

Lemma 2 Let [u, B\ & V he a node for which MDFS performs the operation 
PUSH{[u, B]). Furthermore, at the moment when POP{[u, B]) is performed, only 
strongly simple paths have been constructed by MDFS. Let [x. A] E V' be a node 
such that at the moment when PUSH{[u, B]) is performed, there is a strongly simple 
path P = [u, B],[v, A\,Q,[x, A] with [^,X],[^,X] ^ K for all [z,X] G P. Then 
PUSH{[x, A]) has been performed before POP{[u, B]). 

Remark: Lemmaimplies that either PUSH([x,A]) and POP([x,A]) have been 
performed before the execution of PUSH([n, R]), or both operations have been 
performed between the operations PUSH([n,i?]) and POP([n, R]). 

Proof: Assume that the assertion does not hold. We consider a path P = 
[u, B],[v, A],[v', B],Q',[x, A] of shortest length such that PUSH([x,A]) has not 
been performed before POP([m, B]). It is clear that the edge e = {[u, B], [n. A]) has 
been considered before the execution of POP([M,i?]). By assumption, [n. A] ^ K 
and [n, B] ^ K the moment when e is considered. Hence, MDFS is in Case 2.3. 

If the operation PUSH([n,A]) is performed according to this consideration of 
edge e, then PUSH([n', R]) would be the next operation performed by the algo¬ 
rithm. By assumption, P is a shortest path such that the assertion is not fulhlled. 
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Therefore, PUSH([x, A]) has been performed before POP([n', B]), and hence, before 
POP([m,-B]), a contradiction. 

If the operation PUSH([n, A]) is not performed according to this considera¬ 
tion of edge e, then MDFS is in Case 2.3.i and performs the corresponding while- 
statement. Consider the moment when MDFS terminates this while-statement; 
i.e., T[,;,a] = 0- Let [z,A] be the hrst node on P for which PUSH([z,H]) has 
not been performed. Since [x. A] has this property, the node [z, A] exists. Let 
P = [m, H], [n, H], Qi, [?/, H], [z, H], Q 2 , [x, A]. By construction, each node on Pi = 
[u, B], [n. A], Qi, [y, B] is pushed and each edge on Pi is considered. Hence, all these 
nodes and edges are in Texp such that Pi is formed by MDFS. Furthermore, the 
edge {[y,B], [ 2 , H]) is a considered weak back edge. Therefore, MDFS has found 
the path Pi, [z, A]. By assumption, PUSH([z, H]) has never been performed. Since 
[z,B] ^ K when PUSH([ti, P]) is performed it holds that POP([ 2 ;,P]) has been 
performed. Hence, [z,A\ G and hence, L]^v,a] 7^ 0- But this contradicts 

L[v,a] = 0 such that the lemma is proved. ■ 

The following lemma shows that in a certain situation some nodes are already 
pushed. For the proof of the lemma, we need the notation of the so-called back-path 
r(S') of a path S in Gm- Essentially, the back-path of S is obtained by changing the 
direction of the edges on S and running S backwards. To get a formal dehnition 
for r(S'), we denote for w G V' 




t a w = s 

s if w = t. 


Let S = Wi, W2, ..., tCfc be a path in Gm- The hack-path r{S) of S is dehned by 


r{S) := r{wk),r{wk-i), - - • ,r(wi). 


The following lemma is a direct consequence of Lemma 

Lemma 3 Let [m, B] E V be a node for which MDFS performs the operation 
PUSH{[u, B]). Furthermore, at the moment when POP{[u, B]) is performed, only 
strongly simple paths have been constructed by MDFS. If there exists a strongly sim¬ 
ple path P = [v,A\,Q, [w,B] such that at the moment when PUSH{[u, B]) is per- 



Em then for all [z,X] G P, the operations PUSH{\z,X]) and PUSH{[z,X]) have 
been performed before the execution of the operation POP{[u, B]). 

Proof: For [z, X] G P, the assertion follows by an application of Lemma 2 to the 
path [u,B],P. With respect to [z,X] where [z,X] G P, we apply Lemma 2 to the 
path [u, P], r(P). ■ 
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Note that by the dehnition of |-^[•u,A]| > 0 implies that PUSH([m, A]) 

and POP([m, A]) have been performed. The following lemma extracts properties of 
the algorithm MDFS which enable us to prove the correctness and to develop an 
efficient implementation of MDFS. 

Lemma 4 MDFS maintains the following invariants: 

1. MDFS constructs only strongly simple paths. 

2. < 1? for all [w,A\ G V . 

3. Assume that the algorithm performs the assignment := Then 

after the execution of PUSH{[u, A]) always 

Remark: Invariant 2 and Invariant 3 are not needed for the correctness proof. 
But we shall need these invariants for an efficient implementation of the algorithm. 
Moreover, the proof of Invariant 1 is easier if we prove all invariants simultaneously. 

Proof: Consider the hrst situation in which one of the three invariants is not 
maintained. Three cases have to be considered. 

Case 1: Invariant 1 is not maintained. 

Only a PUSH-operation can destroy the property “strongly simple”. Note that 
a PUSH-operation cannot affect Invariant 2 or Invariant 3. Lemma 1 implies that 
this PUSH-operation occurs during the consideration of an edge e = {[v,B], [w, H]). 

If [w,A\ is not marked “pushed”, then Case 2.3.ii applies and PUSH([t(;,H]) 
is performed. Since [w, B] ^ K, the only possible situation in which this PUSH- 
operation destroys the property “strongly simple” is the following: On a current 
search path there is a subpath Q which is caused by an application of Case 2.3.i 
of MDFS such that [w, B] G Q. Hence, there exists [u, A] G V such that the 
addition of Q to this current search path is caused by the operation PUSH([m, H]). 
By construction, the assumptions of Lemma are fulfilled with respect to [u, B] 
and [w, H] on P. Hence, by Lemma PUSH([t(;, H]) has been performed before 
POP([m,R]), and hence, before PUSH([m,H]), a contradiction. 

Hence, [ta. A] is marked “pushed” such that Case 2.3.i of MDFS applies. By 
Invariant 2, |L[iu,a]| < 1- We thus write = [w, instead of 

By Case 2.3.i, for the node [u, A] = the operation PUSH([m, H]) is performed 

such that MDFS extends the current search paths by a path [tc,H],(5, [m,H], but 
only [u, A] is pushed. By the definition of and by Lemma 3, the opera¬ 

tions PUSH([z,X]), POP([z,X]), PUSH([z,X]),’and POP([z,X]) have been per¬ 
formed, for all [z,X] ^ Q such that none of these nodes is contained in the current 
MDFS-stack K. Hence, the only possible situation in which PUSH([m,H]) de¬ 
stroys the property “strongly simple” is the following: There is a node [p, X] G 
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[w, A], Q, [u, A], and a subpath Q' of a current search path which is caused by an 
application of Case 2.3.i such that \p,X] G Q' or [p,X] G Q'. Since one end node 
of an edge in the current matching uniquely determines the other end node, we can 
choose [PiX] such that [p, A] G Q'. 

Consider the node [u\A\ G K such that PUSH([m', A]) is the operation which 
adds the subpath Q' to this current search path. Therefore, immediately before 
the execution of PUSH([m', A]), T[p,A] = Hence, by Invariant 3, after the 

execution of PUSH([m', H]), always L^,a] = L\^u',a]- By the choice of [p, H], = 

and hence, = [u,A\ in the situation under consideration. Hence, 

POP([m',H]) has been performed such that H] ^ K, a contradiction. 

Case 2: Invariant 2 is not maintained. 

Then there exists [tc,H], [pi,H], [p 2 ,H] G V' such that = {[Pn^]} before 

the execution of POP([p 2 ,H]) and Li^^a] = {[pi, H], [p 2 , H]} after the execution 
of POP([p2,H]). Hence, MDFS has found a path Pi = [pi, B],Q,[pi, A] with 
[w, A] E Q and found a path P 2 = [p 2 , B], Q', [p 2 , A] with [w, A] G Q'. 

If MDFS has found the path P 2 after the execution of POP([pi, H]), then [tc. A] 
can only be added to Q' in the following way: An operation PUSH([m, A]), caused 
by an application of Case 2.3.i with respect to a node [n. A] (i.e., [m. A] G is 

performed such that the current search path is extended by a path [n. A], Q, [u, A] 
with [tc. A] G Q. But then, [m. A] G T[tu,A] before the execution of PUSH([m, A]). 
PUSH([m, A]) is performed after POP([pi,H]). Hence, [m,A],[pi,A] G be¬ 

tween the execution of these two operations. This contradicts the assumption that 
we consider the situation in which Invariant 2 is not maintained for the first time. 

Hence, MDFS has found the path P 2 before the execution of POP([pi,H]). 
Note that \pi,B] ^ Q'. Otherwise, by Lemma 3, PUSH([pi, A]) is performed before 
POP([p2,H]), and hence, [pi. A] ^ L[w,a] after POP([p 2 ,H]). Let [r. A] be the first 
node on Q' such that [r. A] G Q or [r,B] G Q. Since node [tc. A] has this property, 
the node [r. A] exists. Let 


Q' = Qi: [a A]: Q2 a’ld Q 


Qi,[r,A],Q 2 ii[r,A]EQ 
Qi, [r,B],Q 2 if [r,B]EQ. 


Consider the path 


/ Q'l, [r, A],Q 2 , [pi. A] if [r. A] E Q 

1 Q'l, V,A,^{Qi), if V,B] E Q. 


Then Lemma 2 applies with respect to [p 2 , B], [pi. A] and the strongly simple path 
R. Hence, PUSH([pi,A]) is performed before POP([p 2 ,H]), and hence, [pi,A] ^ 
after POP([p 2 ,i?]), a contradiction. 

Case 3: Invariant 3 is not maintained. 
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After the execution of PUSH([u, A]), it holds that L]^w,a] = L\u,a] = 0- We shall 
prove that L]^w^a] = L]^u,a] after the next POP-operation which changes or 

L[u,a]- Then, the assertion follows because of Invariant 2 and the transitivity of 
the relation =. Let POP([p, S]) be the next POP-operation which enlarges Lyw,A] 
or L[n,A]- K]^w,a] denotes the current MDFS-stack, directly after the execution of 
PUSH([r(;, A]). Let K' = Note that [u,B] G \ K'. According 

to the location of [p,-B] with respect to and to K]^u,a]) we have to discuss 

three cases. 

By construction, [p, S] ^ \ K'. Otherwise, POP([p, S]) would be per¬ 

formed before PUSH([n, A]). 

Assume that [p, B] G K]^u,a] \ K' ■ Let [g, B] be the hrst node in iP[u,,A] \ K' such 
that [g, A] G K]^u,a] \ Node [g, S] exists since [u,B] has the property that 

[u,B] G \ K' . Consider the back-path of the path from node [p, S] to node 

[g. A], This back-path implies that [g, B] and [p. A] fulhll the assumptions of Lemma 
2. Hence, PUSH([p, A]) occurs before POP([g, H]). Since [g, 5] G \K\ the 

operation PUSH([p, A]) is also performed before POP([p, H]). Hence, POP([p, H]) 
can enlarge neither L[^,a] nor L[„,a]- 

It remains to consider [p, B] G K'. Let [g, B] G K' be the node nearest to 
the top of K' for which PUSH([g, A]) has not been performed at the moment 
when MDFS performs PUSH([m,A]). Since [p, H] has this property, [g, H] exists. 
By consideration of the back-path of the path from [g, H] to it is easy to 

prove that MDFS hnds a path from [n. A] to [g. A] not containing [g,H]. Hence, 
L[u,a] = b) after the execution of POP([g, H]), and hence, [q^B] = [p, H]. Since 
MDFS has found a path from [tc. A] to [u. A] which does not contain [g, H], it holds 
that L[^,a] = [g,^] = [p, A]. ■ 

Now, the correctness of the algorithm MDFS can easily be derived from Lemma 
|2]and Lemma IH 


Theorem 3 MDFS constructs a strongly simple path from s to t iff such a path 
exists. 


Proof: Assume P = s, [uq, H], [ui. A], [n'^, H],..., H], [n^, A], f is a strongly 

simple path from s to t. It is clear that MDFS considers the edge (s, [ug, B]) and per¬ 
forms the operation PUSH([nQ, B]). (Note that Vq is M-free.) Hence, [ug, B], [n^. A] 
fulhll the assumptions of Lemmaj^with respect to the path [ug, B], [ui. A],..., [n(._i, B], 
[vr,A]. By Lemma[^ MDFS performs PUSH([nr, A]) and hence, PUSH(f). There¬ 
fore, MDFS constructs a path from s to t. By Invariant 1 of Lemma MDFS 
constructs only strongly simple paths. | 
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2.5 An implementation of MDFS 

Now we shall describe how to get an efficient implementation of MDFS. Only two 
parts of the algorithm are nontrivial to implement. 

1. The manipulation of L]^w,a]) [wiA\ G V. 

2. The reconstruction of a strongly simple path P from s to t which is con¬ 
structed by the algorithm. 

For the solution of both subproblems it is useful to perform the POP-operations not 
explicitly and to maintain the whole MDFS-tree T. This can be done as follows: 
The data structure is a tree T. A pointer TOP always points to TOP(iF) in T. 
The current MDFS-stack K is represented by the unique path from the root s of T 
to TOP(iF) in T. For performing the operation POP, the pointer TOP is changed 
such that it points to the unique direct predecessor in T. When we perform a 
PUSH-operation, the node in T to which TOP(iF) points obtains a new leaf. After 
the PUSH-operation TOP(iF) points to this new leaf. 

Invariant 2 and Invariant 3 are the key for the efficient implementation of MDFS. 
Now we shall describe the update of Li^ui,a]- By the dehnition of we only 

have to change after a PUSH- or after a POP-operation. More exactly, we 

have to perform after PUSH([m, A]) the operation := 0 if = ["W, A] and 

after POP([M,i?]) the operation := [w. A] if PUSH([m,A]) has never been 

performed and MDFS has found a path P = [tc, A],Q, [u, A] for which [u, B] ^ Q. 

After the execution of POP([m, B]), if PUSH([m, A]) has never been performed, 
MDFS needs all nodes [tc. A] for which a path P = [tc. A], Q, [m. A] such that 
[u, B] ^ Q has been found by MDFS. This can easily be done by any graph search 
method like depth-first search on the current starting at node [u, A] and run¬ 
ning the considered edges backwards. When the node [u, B] is reached, a backtrack 
is performed. But with respect to efficiency, it is useful to investigate the properties 
of MDFS and to rehne the backward graph search. 

First, we shall characterize the paths P = [tc. A], Q, [m. A] with [u,B] ^ Q, 
found by MDFS. Let P = ci, 62 ,..., e^. Then, the following properties are fulhlled: 

1 . Ct is a weak back edge. 

2. If we start in edge e* and consider P backwards, then we see a nonempty se¬ 
quence of tree edges followed by a single cross, forward or back edge, followed 
by a nonempty sequence of tree edges followed by a single cross, forward or 
back edge, and so on. 

Hence, after the execution of POP([m, P]), we need the following sets of nodes: 
R[u,a] ■= {['y, B] E V' \ ([n, P], [u, A]) is a weak back edge} 
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and for some [q-,^ G V 

E[q,A] ^ I b)^]) is a cross, forward, or back edge}. 

According to Invariant 3, dnring the backward search, some snbpaths can be 
skipped over. Therefore, we need the following set of nodes 

D[qA] ■= e I ^bA] = previonsly}. 

By Invariant 3, T*[q,A] ^ -D[g',A] implies L\^q^A] = -^[g'A]- We say that T*[q,A] is current 
if T*[g,A] ^ ^lq\A], for all [q',A] E V' \ {[g, A]}. According to Invariant 3, we can 
compnte L^^a] in the following way. 

1 . Compnte [q,A] snch that [p. A] G D[q,A], and T*[g,A] is cnrrent. 

2. If [q, A] does not exist, then L^^a] = 0- Otherwise, 

r _ / if PUSH([g, A]) has never been performed 

} 0 otherwise. 

As described above, a correct manipulation of the current sets T*[q,A] allows the 
solution of the hrst subproblem. Note that by Invariant 2 of Lemma each 
[p. A] G V' is contained in at most one current set T*[g,A]- If during the backward 
search a node [p. A] is met for which L^^a] 7^ 0 previously, some nodes can be 
skipped over. Hence, we have also to know if L^p^A] 7 ^ 0 previously. This will be 
realized by the correct update of the following set 

L := {\p, A] e V \ Lip^a] 7^ 0 previously}. 

Now we can give a detailed description of the backward search which will be per¬ 
formed after POP([M,i?]). The consideration of those paths P = [w, A],(5, [n. A] 
with [u, B] ^ Q is done in several rounds. In the hrst round, we construct back¬ 
wards all paths without any cross, forward, or back edge. In the second round, 
all paths with exactly one such edge are constructed implicitly, and so on. Let 
denote the nodes in Tf.xp which have been already considered during the backward 
search. If a node in is considered again, the search has not to be continued 
at that node. In the hrst round, we consider the weak back edges {[v,B], [u. A]). 
In the All round, i > 1, we consider those edges {[v,B], [g. A]) G i?[q,A] for which 
L[g^A] = [u, A] is computed in the (i — l)st round. Starting in node [n, B], we follow 
backwards the tree edges as long as a node in U {[u, B]} is reached. If we reach 
a node [p. A] G L, then we compute the current T*[r,A] such that [p. A] G T*[r,A] and 
we jump to [r, B] for the continuation of the backward search. Since we perform a 
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backward search, [r, B] is switched to [r, A\. According to Invariant 3, L\x,a] = 
and hence, L]^x,a] = [«, A] for all [x. A] G -D[r,A]- 

For the organization of the backward search, we use a queue Q which contains 
the start nodes of the next round. During Round i, the start nodes of Round i + 1 
are added to the end of the queue. Since the data structure is a queue, all start 
nodes of a round are treated before the hrst start node of the next round is taken 
away from Q. 

For the reconstruction of a strongly simple path from s to f constructed by the 
algorithm, we have to know the non-tree edges used on the path. Hence, for all 
[r. A] G V we use a variable P[r,A] to store the needed information with respect to 
the node [r. A]. This means, we store in variable P[r,A] that non-tree edge which 
concludes the block of tree edges which contains the tree edge with end node [r. A] 
at the moment when L\^x,a] 7^ 0 for the first time during the backward search. 

The implementation of MDFS must be done with attention to the correct ma¬ 
nipulation of the sets R[q,A]i and T^[ 5 ,a]- The following table describes in 

terms of the case of MDFS, and in terms of the operation which is performed, how 
MDFS has to update these sets. 


ease, operation 

set updating 

Case 1 

no update 

Case 2.1 

• PlwjA] C 

Case 2.2.i 

• PlwjA] C 

Case 2.2.ii 

Case 2.3.i 

R[w,A] • R[w,A] C 

L[w,a] 7^ 0 

no update 

L[w,a] = 0 

:= E[w,a] U {[n, B]} if [w, A] ^ L 

Case 2.3.ii 

no update 

PUSH([m,A]) 

no update 

POP([n,R]) 

EIv,a] ■= {[P; A] MDFS has found a path from 
[p. A] to [n. A] not containing [v,B]} 


In Case 2.1, it is clear that [tc,A] ^ L since POP([tc, A]) has not been performed. 
In Case 2.2.i, [tc. A] ^ L follows directly from [w,B] G K and Lemma 1. Note 
that in Case 2.3.i, subcase T[u,,a] 7 ^ 0, we have to store the information that edge 
{[v,B], [tc. A]) is used. In the implementation, we accomplish this by adding the 
edge ([n, R], [tc. A]) to the node [v,B] in K. Then we obtain the expanded node 
{{[v,B], [tc. A]); [v,B]). The considerations above lead to the following implemen¬ 
tation of the procedure SEARCH. 


procedure SEARCH; 

if top(a:) = t 


18 





then 

reconstruct a strongly simple path P from s to t 
which has been constructed by the algorithm 

else 

mark TOP (if) “pushed”; 

for all nodes [w, Y] e N[TOP{K)] 

do 

iiY = B 

(Case 1) then 

FVSR{[w,B]); 

SEARCH 

(Case 2) else 

if [w, A] G if 

(Case 2.1) then 

E[w,A] ■= 0 {TOP(if)} 

else 

if [w, B]eK 


(Case 2.2) 

then 


if [w,A] is marked “pushed” 

(Case 2.2.i) 

then 


E[w,A] ■= 0 {TOP(if)} 

(Case 2.2.ii) 

else 


R[w,A] ■= R[w,A] U {TOP(if)} 


fi 

(Case 2.3) 

else 


if [tc,A] is marked “pushed” 

(Case 2.3.i) 

then 


if -h[u>,A] 7^ 0 

then 

expand TOP(if) in if to 
((TOP(if),[n;,A]);TOP(if)); 
PUSH(L[^,^]); := 0; 

SEARCH 

else 

if [w, A] ^ L 

then 

E[w^a] ■= E[^,a] U {TOP(if)} 

fi 

fi 

(Case 2.3.ii) else 

PUSH([n;,A]); 


19 



SEARCH 


fi 

fi 

fi 

fi 

od; 

(* let T0P(A:) = [v,B]*) 

Lcur • 

-Dw •= 0; 

Q;= 0; 

for all [g, B] G R[v,a\ 

do 

CONSTRL(([g,R],[t;,A]),[t;,R]); 

od; 

while Q ^ 0 
do 

remove the front node [/c, A\ from Q; 
for all [g, B] G which have not already been 

considered during the backward search 

do 

CONSTRL((|g,B],|i,^]),[t,,B]) 

od 

od; 

POP; 

POP 


CONSTRL is a call of the following procedure. The variable Pcur contains always 
the non-tree edge which concludes the current block of tree edges. 

procedure CONSTRL(([g, R], [m, A]), [x, R]); 

Pcur ■— ([o')-S], [it, A]), 

[z,B] := [g,R]; 

while no node in U {[x, R]} is reached 

do 

starting in the node [z, R], perform a backward search until 
a node in L U U {[x, R]} is reached on the tree edges; 
if [y, A] ^ L is met during the backward search 

then 

;= Rl_ U{[|/,A]}; 

L:=LU{[y,A]}- 
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P[y,A\ ■ Pcuri 

add the node [y, A\ at the end of Q 

fi; 

if [y, A] E L is met by the backward search 

then 

(* Let -D[r,yi] be the current set containing [y, A], *) 
Dl •= -Dl a D\^ 

\z,B] := [r,B] 

fl 

od. 


The reconstruction of a strongly simple path P from s to t constructed by the 
algorithm remains to be explained. Beginning at the end of P, such a path P 
can be reconstructed by traversing the MDFS-tree T backwards. Note that TOP 
points to the end of P, and that the father of each node in T is always unique. 
As long as we traverse tree edges of the algorithm MDFS, we have no difficulty. 
But every time when we meet a node [u, A] which has been added to P by an 
application of Case 2.3.i, we have to reconstruct a subpath [w,A\,Q, [u,A\ which 
has been joined to P. In this situation, the considered portion of T is the expanded 
node (([n, P], [tc, A]); [n, P]); i.e., the structure of T tells us that MDFS has applied 
Case 2.3.i. It remains the reconstruction of the subpath [tc, A], Q, [m, A]. For doing 
this, we start in the node [tc,A]. We use P[w,a] to obtain the non-tree edge of 
MDFS, which finishes the block containing the tree edge with end node [w, A]. Let 
Plw,A] = [y,A]). Then denotes [x,B], and P[^^a] denotes [y,A]. First 

we reconstruct the block from the node [tc. A] to the node = [x,B]. Then 

we reconstruct the block from the node Pj^ = [y, A] to the node P^^ ^j, and so on 
until the node [u, A] is met. Each block can be reconstructed as the path P itself. 
These considerations lead to the following procedure for the reconstruction of an 
augmenting path, constructed by the algorithm. 

procedure RECONSTRPATH(t, s); 

NODEcur ■= t; 

while NODEcur 7^ s 

do 

if father (NODEcur) is not expanded 

then 

NODEcur '■= father (NODEcur) 
else (* let father (NODEcur) = {([v,B], [tc. A]); [v,B])*) 

RECOmTRQ(NODEcur, [w, A]); 

NODEcur ■= b,P] 


21 



fi 


od. 


RECONSTRQ is a call of the following procedure. 

procedure RECONSTRQ([«, A], [w, A]); 
ST :=[w,A\, 

RECONSTRPATH(Fj3., ST); 

while Pgrp ^ [u,A\ 

do 

ST := 

RECONSTRPATH(P^y, ST) 

od. 


The correctness of the manipulation of L]^w,a]-i and the correctness of 

the reconstruction of the M-augmenting path P follow from Lemma and are 
straightforward to prove. The procedure RECONSTRPATH resembles standard 
recursive methods used for the reconstruction of augmenting paths (see e.g. H). 
The time and space complexity of our implementation of MDFS remain to be 
considered. It is easy to see that the time used by the algorithm MDFS is bounded 
by 0{n + m) plus the total time needed for the manipulation of the sets P[g,A]; 
[g, A] G V. If we use linear lists for the realization of the sets P[q,A] with a pointer 
to the node [g, A] for each element of P[g,A ]5 the execution time for each union 
operation is bounded by 0{n). Following the pointer corresponding to [p, A], we 
can End the set containing [p, A] in constant time. At most n union operations are 
performed by MDFS. Hence the total time used for the manipulation of the sets 
D[q^A] is bounded by The time needed for the n union operations can be 

reduced to O(nlogn) if we use the following standard trick, the so-called weighted 
union heuristic. 

We store with each set the number of elements of the set. A union operation is 
performed by changing the pointer of the smaller of the two sets which are involved 
and updating the number of elements. Every time when the pointer with respect 
to an element is changed, the size of the set containing this element is at least twice 
of the size of its previous set. Hence, for each element, its pointer is changed at 
most logn-times. Hence, the total time used for all union operations is 0{n\ogn). 
Altogether, the total time used for the augmentation of one augmenting path is 
0{m + n\ogn). 

If we use for the update of the sets D\^q^A] disjoint set union [M], the total 
time can be bounded to be 0{{m + n)a{m,n)) where a is the inverse Acker ma un 
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function. Note that for each node [p, A] one hnd operation suffices for the decision 
of L[p^a]. Furthermore, we can reduce these bounds to 0{m + n) using incremental 
tree set union [12]. The space complexity of MDFS is bounded by 0{m + n). The 
considerations above lead to the following theorem. 

Theorem 4 MDFS can be implemented such that it uses only 0{m + n) time and 
0{m + n) space. 


3 The Hopcroft-Karp approach for general graphs 

In 1973, Hopcroft and Karp [2l] proved the following fact. If one computes in one 
phase a maximal set of shortest pairwise disjoint augmenting paths and augments 
these paths then O(yTf) such phases would suffice. In the bipartite case, they 
have described an elegant simple 0{m + n) implementation of an entire phase. 
Let us sketch this implementation. First they have reduced the problem of hnding 
augmenting paths to a reachability problem in a directed graph Gm with two 
additional nodes s and t. Then, by performing a breadth-first search on Gm with 
start node s until the target node t is reached, they have obtained a layered, directed 
graph Gm for which the paths from s to f correspond exactly to the shortest 
M-augmenting paths in G. Using depth-first search, they find a maximal set of 
pairwise disjoint M-augmenting paths. Whenever an M-augmenting path is found, 
the symmetric difference is applied to the path and the current matching, the path 
and all incident edges are deleted and the depth-hrst search is continued. Breadth- 
hrst search and depth-hrst search take 0(m+n) time. Hence, the implementation of 
Hopcroft and Karp of a phase has time complexity 0{m + n). Since M-augmenting 
paths can be found in general graphs by a slightly modihed depth-hrst search 
(MDFS), the following question suggests itself: Can we get an implementation 
of an entire phase by performing something like breadth-hrst search followed by 
MDFS? We shall give an affirmative answer to this question. 

3.1 The description of a phase 

Let G = (U, E) be an undirected graph, M be a matching of G, and Gm = (U', Em) 
be the directed graph as dehned in Section 2.2. Our goal is to construct from Gm 
a layered directed graph Gm = (U', Em) such that 

1. the Ith layer contains exactly those nodes [u,X] G V such that a shortest 
strongly simple path from s to [u,X] in Gm has length /, and 

2. Gm contains all shortest strongly simple paths from s to f in Gm- 
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The level of a node [n,X] G V is the length of a shortest strongly simple path 
from s to [n,X]. In Gmi the *th layer contains exactly the nodes of level i. It is 
clear that s is the only node in Layer 0. By the structure of Gm, level([n, 5]) is 
odd and level([n, A]) is even for all v & V. Since breadth-hrst search (BPS) on Gm 
with start node s hnds shortest simple distances from s and not shortest strongly 
simple distances, BPS cannot be used directly for the construction of Gm- But we 
can modify BPS such that the modihed breadth-hrst search (MBPS) hnds shortest 
strongly simple distances. Remember that for the construction of the (/-|-l)st layer, 
BPS needs only to consider the nodes in Layer /, and to insert into the (/ -|- l)st 
layer all nodes w which fulhll the following properties: 

1. There is a node v in the Ith layer with {v,w) G E. 

2. Level(t(;) has not been dehned. 

With respect to hnding strongly simple distances from s, the construction of the 
(/ -|- l)st layer is a bit more difficult. By the structure of Gm-, the level of a free 
node [w,B\ is one and the level of a non-free node [w,B] is well-dehned by the 
level of the unique node [v,A\ with ([n, A], [w,B]) G Em- Hence, the construction 
of odd layers is trivial. Por odd /, we shall describe the construction of the (/ -|- l)st 
layer under the assumption that Layers 0,1,2,..., / are constructed. It is clear that 
similar to BPS, MBPS has to insert into the (/ -|- l)st layer all nodes [w, A] G V 
which fulhll the following properties: 

1. There is a node [n, H] in Layer I with {[v,B], [tc,H]) G Em, and there is a 
strongly simple path from s to [n, B] of length / which does not contain the 
node [tc, B]. 

2. Level([t(;, H]) has not been dehned. 

But these are not all nodes which MBPS has to insert into Layer Consider the 
example described by Pigure 2. Note that level([v 7 , R]) = 7 but level([v 3 , H]) 7 ^ 8 , 
since the unique shortest strongly simple path from s to [^ 7 , B] contains [^ 3 , B]. The 
unique strongly simple path P from s to [^ 3 , A] has length 14. Hence, level([n 3 , H]) = 
14. Therefore, MBPS has to insert nodes [w, A] G V into Layer 1 + 1 for which there 
is a shortest strongly simple path P = s,[vi, B],..., [vi, B], [tc. A] with level([n«, B]) < 
1. Por the treatment of these nodes and for the knowledge if there is a strongly 
simple path from s to [n, B] of length I which does not contain the node [w, B], the 
following notation is useful. 

Let T C V such that level([n,X]) is dehned for all [v,X] G T. We denote by 
DOM(T) the set of those nodes [m, B] eV'\T which satisfy: 

a) Por all [v,X] G T, all shortest strongly simple paths from s to [n, W] contain 
[m,R], 
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[Vb, B] -► [ve, A] - [vn, B] 



s 


Figure 2: Further node which has to be inserted into Layer I + 1. 

b) level([n, A]) has not been dehned, and 

c) level([ta,-B]) < level([M,-B]) for all [w,B] G V' satisfying a) and b). 

If such a node [u, B] does not exist then DOM(T) denotes the node s. Furthermore, 
DOM({s}) denotes the node s. We shall use DOM(T) only for subsets T of V' of 
size at most two. Next we shall show that always |DOM(T)| = 1. This will be a 
direct consequence of the following lemma. 

Lemma 5 Let P = s,[vi, B],[v 2 , A\,... ,[vi, Z] he any shortest strongly simple 
path from s to [vi,Z], i.e., levelifvi, Z]) = 1. Let [vj,X] E P be a node with 
level{[vj, X]) > 1. Then level([nj, X]) = j and level([ni, X]) < level([nj, X]) for all 
i < j. 

Proof: Since [t^j,X] is the jth node on the strongly simple path P it holds that 
level([nj,X]) < j. Suppose level([nj,X]) < j and let Q,[vj,X] be any shortest 
strongly simple path from s to [vj,X]. Let P 2 := [uj+i, X], [^^+ 2 , X],..., [n^, Z]. 
Note that all nodes on P 2 are not M-free. Let R := Q, [vj, X], P 2 . By construction, 
|B| < |P|. Since P is a shortest strongly simple path from s to [vi,Z], the paths 
Q and P 2 cannot be strongly disjoint. Let [z,Y] be the hrst node on Q such that 
[z.,Y] or [z,Y] is on P 2 . If [z.,Y] is on P 2 then Y = A. \iY = B then the unique 
node [x,A\ with ([x, A], [z,B]) G Em has to be the direct predecessor of the node 
[ziB] on both paths Q and P 2 . This contradicts the choice of the node [2:,X]. 
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If [z.,Y] is on P2 then also Y = A. \i Y = B then the node [x,A\ such that 
([x, y 4 ], [z, B]) G Em has to be on Q the direct predecessor of [z, B] and [x, B] has 
to be on P2 the direct successor of [z, A]. But this would also contradict the choice 
of the node This shows that in any case Y = A. Let 


Q = Qi, [z, A], Q2 and P2 


P2u[z,A],P 22 ii[z,A]EP2 

P21,[z,B],P 22 if[z,B]eP2. 


If [z,A] G P2, then S := Qi, [z^A\^P22 would be a strongly simple path from s 
to [vi,Z] shorter than P, a contradiction. Hence, [z^B] G P2. But then S : = 
Qi, [z, H], r(P2i), [vj,X] would be a strongly simple path from s to [vj,X] shorter 
than 1 . This contradicts level([nj,X]) > /. Hence, level([nj,X]) = j. 

Since [vi,Y] is the ith node on P it holds that level([nj, X]) < i. Hence, i < j 
implies level([uj, X]) < level([nj, X]). | 

The following lemma is a simple consequence of Lemma 

Lemma 6 Let T C V, T 7^ 0 such that level([n,X]) is defined for all [n,X] G T. 
Then the following statements hold true: 

a) |DOM(T)| = 1 . 

b) Let DOM(T) = [m,P]. Then after the definition 0/level([ m, H]), always 
DOM(T) = DOM([m,P]). 

Proof: a) Assume that |DOM(T)| > 1 . Let [mi,P] and [u2, B] be two dis¬ 
tinct elements of DOM(T). By the consideration of any shortest strongly sim¬ 
ple path from s to [v,X] for any [v,X] G T, applying Lemma we obtain 
level([Mi, P]) < level([M2, P]) or level([M2, P]) < level([Mi, P]). This contradicts 
Part c) of the definition of DOM(T). 

b) After the definition of level([M, A]), level(DOM(T)) decreases. Hence, by Lemma 

DOM(T) has to be below [u,B] on all shortest strongly simple paths from s to 
[n, X] for all [n, X] G T. Hence, DOM([m, P]) is on all shortest strongly simple paths 
from s to [n,X] for all [n,X] G T. Assume that DOM(T) 7^ DOM([m, P]). Then 
level(DOM(T)) > level(DOM([M, P])). By the dehnition of DOM([m,P]), there is 
a shortest strongly simple path Pi, [u,B] from s to [u,B] which does not contain 
DOM(T). Hence, Pi, [u, B] cannot be a subpath of a shortest strongly simple path 
from s to [n,X] for a node [n,X] G T. Let Qi, [u,B],Q2 be any shortest strongly 
simple path from s to [n,X]. Lemma [^implies that |Pi,[m,P]| = \Qi,[u, B]\. Since 
Pi, [u,B],Q2 is not a shortest strongly simple path, the paths Pi and Q2 are not 
strongly disjoint. Let [z. A] be the first node on Pi such that [z, A] or [z, B] is on Q2- 
Exactly as in the proof of Lemmawe prove that such a node cannot exist. Hence, 
level(DOM(T)) < level(DOM([M, P])). Therefore, DOM(T) = DOM([m,P]). | 
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Now we shall describe a phase in detail. Note that to each node w ^ V there 
correspond two different levels, namely level([w, A]) and level([w, i?]). The first 
level of w denoted by leveli(w) is the smaller one of these two levels. The other 
level denoted by level2(tc) is the second level of w. 

In [SI E], we have constructed during the (/ + l)st phase exactly the shortest 
strongly simple paths of length / + 1 . Especially for the treatment of the shortest 
strongly simple paths P = s, [ui, B],..., [vi, B], [w, A] with level([nj, B]) < I, we 
have used some sophisticated data structures to achieve a certain time bound. 
In |S], we have used the dynamic nearest common ancestor algorithm of [IB] in 
combination with Fibonacci heaps na to get an O(mlogn) implementation of a 
phase. We have also described an alternative to the use of Fibonacci heaps which 
uses two-dimensional arrays to get an 0 {m + implementation. A better way 
is to use an idea of Micali and Vazirani [BH IBS]- For all nodes w, they compute the 
first level during Phase leveli(tc) and the second level during Phase ^(leveli(tc) -|- 
level2(w) — 1 ) + 1 . We shall incorporate this idea into our framework. 

Each phase separates into two parts. Both parts use breadth-first search. The 
first level of each node w is computed during Part 1 of Phase leveli(ta). The second 
level of each node w is computed during Part 2 of Phase ^(leveli(tc) -|- level2(tc) — 
1 ) -|- 1 . Before describing the two parts of a phase, we shall investigate the structure 
of shortest strongly simple paths defining the first and the second level of a node. 
First we characterize exactly those shortest strongly simple paths which define the 
first level of a node to be / 1. 

If I is even then leveli(w) = Z -|- 1 for a node w G F iff for the unique node v G V 
such that ([n. A], [w, B]) G Em it holds that level([n. A]) = leveli(n) = 1 . Note that 
level([f, 5 ]) < I implies that level([w. A]) < I such that level([ta, 5 ]) cannot be the 
first level of the node w. Each shortest strongly simple path from s to the node 
[n. A] followed by the edge ([f, A], [w, B]) is a shortest strongly simple path from s 
to [w,B]. Furthermore, there is no other shortest strongly simple path from s to 
[w,B]. 

If I is odd then leveli(w) = Z -|- 1 for a node w G F iff Ievel([ta, A]) and 
Ievel([ta, i?]) are larger than Z and there is a node [v,B] with level^v, i?]) = Z 
and the edge {[v,B], [ta. A]) is in Em- Exactly the shortest strongly simple paths 
from s to such a node [v,B] followed by the edge {[v,B], [w. A]) are the shortest 
strongly simple paths from s to [w. A]. 

Next, we shall characterize exactly those shortest strongly simple paths which 
define the second level of a node to be Z -|- 1 . 

If Z is even then level2(tc) = Z -|- 1 for a node tc G F iff level([tc, S]) > Z, 
level([r(;. A]) < Z, and for the unique node v G V such that ([n. A], [w, 5 ]) G 
Em it holds that level([n,A]) = Z. Note that the node [ta. A] cannot be on a 
strongly simple path from s to [n. A] since the unique direct successor of [w, A] 
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has to be the node [v,B]. With respect to the node [n, A], two cases are possible, 
level([n, A]) = leveli(t;) and level([n, A]) = level 2 (n). Exactly the shortest strongly 
simple paths from s to the node [n, A\ followed by the edge ([n, A], [tc, B]) are the 
shortest strongly simple paths from s to [w, B]. 

If I is odd and level 2 (tc) = / + 1 for a node w E V then level([tc, A]) > I 
and level([t(;, -B]) < /. With respect to a shortest strongly simple path P from s to 
[tc, A], for the direct predecessor [n, B] of [tc, A], exactly two situations are possible, 
level([n,-B]) = I and level([n,-B]) < /. 

In the case that level([v,-B]) = I two subcases can happen, level([n,-B]) = 
leveli(n) and level([n,-B]) = level 2 (n). Exactly the shortest strongly simple paths 
P from s to such a node [n, B] with the property that [tc, -B] is not on P followed 
by the edge ([n, B], [tc, A]) are the shortest strongly simple paths from s to [tc. A]. 

If level([n, -B]) < / then consider any shortest strongly simple path 
P = P', [n, B], [tc. A] from s to [tc. A]. Let [x, B] be the last node on P', [n, B] such 
that the length of the subpath from s to [x,B] of P is equal level([a;, P]). Since 
level([n,P]) < / = level([tc, A]) — 1, the node [x,B] has to be on P'. Since the 
level of the direct successor of s on P is one, the node [x, B] on P' exists. Let 
P' = Pi, [x, P], [y, A], P 2 . As we shall prove later, [y. A], P 2 , [n, P], [tc. A] has to be 
on the back-path of a shortest strongly simple path from s to the node [|/,P]. 

Next, we shall describe MBFS in detail. At the beginning, W = 0 and Em = 0- 
In the first phase, the node s is inserted into Layer 0. For each M-free node w E V, 
the node [w, B] is inserted into Layer 1 and the edge (s, [w, P]) is inserted into Em- 
Assume that I > 0 and Phase / is finished. We shall give a detailed description of 
the two parts of Phase I + 1. 

Part 1 of Phase I + 1 

If I is even then MBFS adds for all nodes [n. A] with level([n. A]) = / = leveli(v) 
the unique edge ([n. A], [w,B]) E Em to Em and inserts the node [w,B] into the 
(/ -|- l)st layer. Note that level([n. A]) = level 2 (n) implies that the first level of node 
w is already dehned. 

If I is odd then MBFS considers all edges {[v,B], [tc. A]) with level([n,P]) = I 
and ([n,P], [tc. A]) has not been considered during Part 2 of a previous phase. If 
the edge ([n, P], [tc. A]) has been considered during Part 2 of a previous phase then 
leveli(t(;) is already dehned. Three cases are possible. 

Case 1: level([t(;. A]) > I and level([tc, P]) > /. 

MBFS inserts the node [tc. A] into the (/-|-l)st layer and adds the edge ([n, P], [tc. A]) 
to Em- 

Case 2: level([tc. A]) > / and level([tc, P]) < 1. 

Then the hrst level of the node w is already dehned. But for the computa¬ 
tion of the second level of w and some other nodes, the edge ([n, P], [tc. A]) is 



needed. Hence, MBPS inserts the pair {[v, B],[w, B]) into the set E{k) where 
k := ^(level([n, B]) + level([tc, B])). 

Case 3: level([t(;, H]) < 1. 

Then the hrst level of the node w is also dehned. But the edge {[v,B], [tc,H]) 
might be necessary for the computation of the second level of some nodes. Hence, 
if level([tc, H]) is already computed then the pair {[v, B],[w, B]) is inserted into 
E{k) where k := ^(level([n, B]) + level([t(;, B]). Otherwise, this pair is inserted into 
E{k) directly after the computation of level([tc, H]). 

These are all cases. Next we shall describe the second part of Phase / + 1. 

Part 2 of Phase I + 1 

During Part 2 of Phase / + 1, level 2 (tc) will be computed for all nodes w with 
|(leveli(tc) + level 2 (tc) — 1) = /. MBPS considers all pairs of nodes [x,Z], [u^E] 
where Z = A if I is even and Z = B if I is odd such that 

1 . level([x, Z]) and level([|/, Z]) have been dehned, 

2. ([x, Z], [y,Z]) e Em, and 

3. |(level([a;, Z]) + level([|/, Z])) = 1. 

Note that these are exactly the pairs in E{1). Starting in [x, Z] and [y, Z], MBPS 
performs a breadth-hrst search on the back-paths of the current layered network un¬ 
til DOM({[a;, Z], [y, Z]}) is reached. All visited nodes [u,X] such that level([M,X]) 
has not been dehned and [u,X] ^ DOM({[a;, Z], [y, Z]}) are inserted into Layer k 
where k = level([a;, Z]) -|- level([?/, Z]) -|- 1 — level([M, X]). The edges {[x,Z], [y.,Z]) 
and ([y, 2’], [x,Z]) are inserted into Em- Purthermore, MBPS adds the traversed 
edges which are not in Em to Em- 

3.2 The correctness proof of MBFS 

We say that a path P is constructed by MBPS if all edges on P are inserted into 
Em- We have to prove the following: 

1. Por all nodes w E V, MBPS computes both levels leveli(t(;) and level 2 (tc) 
correctly. 

2. The layered network Gm = (H', Em) computed by MBPS contains all shortest 
strongly simple paths from s to t. 

The correctness of MBPS is a direct consequence of the following theorem. 
Theorem 5 Let w E V and leveli(t(;) = level([tc,X]). Then the following holds: 
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a) MBPS defines \eve\{[w,X]) eorreetly during Part 1 of Phase \eve[{[w,X]). All 
shortest strongly simple paths from s to [w, X] have been constructed after the 
termination of Part 1 of Phase level([w, X]). 

b) Level{[w, X]) has been defined and all shortest strongly simple paths from 
s to [w,X] have been constructed after the termination of Part 2 of Phase 
^(level([w, A]) + level([w, B]) — 1) + 1. 

Proof: We prove the theorem by induction on the number of performed phases. 
By construction, it is clear that the assertion of the theorem is fulhlled after the 
termination of Phase 1. Note that no second level of a node has to be computed 
during Phase 1. Assume that / > 0 and that the assertion of the theorem is fulhlled 
after the termination of Phase 1. We shall prove that the assertion also holds after 
the termination of Phase / + 1. 

First we shall prove that for all nodes with hrst level / + 1, their hrst level is 
computed correctly during Part 1 of Phase I + 1. Furthermore, the corresponding 
shortest strongly simple paths are constructed after termination of Part 1 of Phase 
/ + 1. No level of any other node is determined during Part 1 of Phase / + 1. 

If the assertion of the theorem does not hold for a node [u, B] then the assertion 
also does not hold for the unique node [n. A] with ([n. A], [u,B]) G Em- Hence, if 
/ is even, the induction hypothesis implies that the assertion is fulhlled after the 
termination of Part 1 of Phase / + 1. Therefore, we have only to consider the case 
that I is odd. 

Consider any node u gV with leveli(M) = level([M, A]) = / + 1. Let P be any 
shortest strongly simple path from s to [m. A]. Let [x,B] be the direct predecessor 
of [m. A] on P. Since level([M, H]) > level([M, A]), the node [u,B] cannot be a 
node on a shortest strongly simple path from s to [x,B]. Hence, level([a;, H]) = 
level([M, A]) — 1. Since the assertion of the lemma has been maintained in previous 
phases, the node [x,B] and the edge {[x, B],[u, A]) have been considered during 
Part 1 of Phase I + 1. Therefore, the node [m. A] is inserted into Layer I + 1 and 
the edge {[x,B], [m. A]) is inserted into Em- Since each shortest strongly simple 
path from s to [x, B] is contained in the current layered network, the path P is 
constructed as well. On the other hand, if we have a shortest strongly simple path 
P from s to a node [x,B] of length I and an edge {[x,B], [m. A]), where leveli(M) 
is not dehned, then the path P, [u, A] is a shortest strongly simple path from s to 
[m. A]. Hence, leveli(M) = 1 + 1. Therefore, only for nodes with hrst level / + 1, the 
hrst level is computed during Part 1 of Phase / + 1. This shows that the assertion 
is fulhlled after the termination of Part 1 of Phase / + 1 such that Part a) of the 
theorem is proved. 

Next we shall show that the assertion of the theorem is fulhlled after the ter¬ 
mination of Part 2 of Phase I + 1. It is clear by construction that during Part 2 of 
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a phase only the second level of some nodes is computed. We have to prove that 
after the termination of Part 2 of Phase / + 1, the following properties are fulhlled. 

1. For all nodes [u,X] such that level([M, X]) := k is computed during Part 2 of 
Phase Z + 1 it holds that |(level([M, A]) + level([M, B]) — l) = l and level([M, X]) 
is computed correctly. 

2. For all nodes u E V with ^(level([M, A]) + level([M, 5]) — 1) = /, the second 
level of u is computed correctly during Part 2 of Phase I + 1 and all shortest 
strongly simple paths from s to [u, X] with level([M, X]) = level 2 (M) have been 
constructed after the termination of Part 2 of Phase / + 1. 

First, we shall prove that the hrst property is fulhlled and then we shall consider 
the second property. 

Assume that the hrst property is not fulhlled. Consider the hrst situation 
such that the hrst property is not fulhlled. Assume that this situation occurs 
during the backward search with respect to the pair [x, Z] [y, Z] of nodes because 
of the assignment level([M,X]) := k. Since a breadth-hrst search on the back- 
paths is performed and this is the hrst time that the hrst property is not fulhlled, 
level([M, X]) has to be larger than k. Furthermore, for each predecessor of [u, X] on 
the back-path, its second level is computed correctly. By the construction of the 
algorithm MBFS, there also holds level([x, Z]) -|- level([j/, Z]) = 21. First we shall 
prove that |(level([«. A]) -|- level([M, 5]) — 1) = I and then that there is a shortest 
strongly simple path from s to [m, X] of length k. 

By construction, k = level([x, Z]) + level([|/, Z]) -|- 1 — level([M, X]). Hence, 

level([M, X]) -|- level([M, X]) = level([a;, Z]) + level([j/, Z]) -|- 1 = 2/ -|- 1 
and therefore, 

-(level([M, A]) -|-level([«, H]) — 1) = /. 

Assume that no shortest strongly simple path from s to [m, X] of length k ex¬ 
ists. By Lemma on all shortest strongly simple paths which contain the node 
[m,X], the length of the subpath from s to [m,X] is equal level([M,X]. Fur¬ 
thermore, since level([M,X]) is computed during the backward search, [m,X] ^ 
DOM({[x,Z],[?/,Z]}). Moreover, level(DOM({[a:, Z], [j/, Z]})) < level([M, X]). 

Since no shortest strongly simple path from s to [m, X] of length k exists, the 
node [m, X] has to be on the all shortest strongly simple paths from s to the direct 
predecessor of [m, X] during the backward search. Each such strongly simple path 
consists of a shortest strongly path from s to [x, Z] or [?/, Z] and the backpath 
of a shortest strongly simple path from s to [y,Z] and [a:, Z], respectively until 
the node ['U,X] is reached without the node [m,X]. Therefore, the subpath from 
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s to [x,Z] and [y,Z], respectively of these paths contains the node [u,X]. Since 
[u,X] 7 ^ DOM({[a;, Z], [y, Z]}), there is a shortest strongly simple path from s 
to at least one of the nodes in {[x, Z], [y,Z]} which does not contain the node 
[m,X]. W.ho.g., let R,[x,Z] be such a path. Let S = S'!, [m, X], S' 2 , [?/, 2’] and 
S' = [m, X], S '21 [x, Z] be any shortest strongly simple paths from s to [y, Z] and 

to [x, Z], respectively such that the path Si, [u, X], S' 2 , [y, Z], [x, Z],r{S' 2 ) is strongly 
simple and constructed during the backward search. Consider the path 


T:= R,[x,Z],[y,Z],r{S,),[u,X]. 


By construction, |T| = k. Since level([M, X]) > k, the paths R and r(S' 2 ) and hence, 
the paths R and S' 2 , are not strongly disjoint. Let [c, A\ be the hrst node on R such 
that one of the nodes [c, A\ or [c, B] is on S' 2 . Let 


R = Ri, [c, A\, i ?2 and S '2 


S21,[c,A],S22 if[c,kl]e^2 

S' 21 , [c, B], S '22 if [c, 5] G S' 2 . 


If [c, B] is on S '2 then i?i, [c, A], r(S' 2 i), [u, X] would be a strongly simple path from 
s to [u, X] shorter than level([M, X]), a contradiction. Hence, [c, A\ is a node on S' 2 . 
By construction, i?i, [c, A\, S' 22 , b, Z] is strongly simple. Hence, 


Ri,[c,A]\ > |,Si,b,X],,S 2 i,[c,H]|. 


If this inequality is strict then [S'!, [u, X], S' 21 , [c, H], i? 2 , [x, Z]\ < level([a;, Z]). Hence, 
the paths S'!, b,X],S' 2 i and R 2 are not strongly disjoint. Let [d,A\ be the hrst 
node on S'!, [u, X], S '21 such that one of the nodes [d, A\ or [d, B] is on i? 2 - If [d, A\ 
is on i ?2 then a strongly simple path from s to [x,Z] shorter than level([a;, Z]) 
can easily be constructed, a contradiction. Hence, [d,B] is a node on i? 2 - Let 
R 2 = R 21 , [d, B], i? 22 - We distinguish two cases. 

Case 1: [d,A] is on S' 21 . 

Let S '21 = S' 211 , [d, A], S' 212 - Consider the path 


T := Ri, [c, A],R 2 i, [d, H],r(S' 2 ii), b,X]. 


By construction, T is a path from s to [u, X] shorter than level(b, X]. By the choice 
of the nodes [c, H] and [d, H], the path T is also strongly simple, a contradiction. 
Hence, Case 1 cannot happen. 

Case 2: [d,A\ is on S'!. 

Let S'! = S'!!, [d, H], S' 12 . Consider the path 


T := 5ii, [d, H], r(i? 2 i), [c, B],r{S 2 i), [u, X]. 
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Note that |T| < level([M,X]). Hence, T cannot be strongly simple. By the choice 
of the nodes [c, H] and [d, H], the paths R21 and S'21 are not strongly disjoint. Let 
[h, A\ be the hrst node on i?2i such that one of the nodes [h, A\ or [h, B] is a node 
on S'21. If [h, A\ is on S'21 then a strongly simple path from s to [u, X] shorter than 
level([M, X]) can easily be constructed, a contradiction. Hence, [h,B] is a node on 
S'21. Let i?2i = d?2ii, [h, H], i?2i2 and S'21 = S'211, [h, H], S'2i2- Consider the path 

T := /?!, [c, H], i?2ii, [K dl], r(^2ii), [w, X]. 

By construction, T is a path from s to [m,X] shorter than level([M,X]). By the 
choice of the nodes [c, H] and [h, H], the path T is also strongly simple, a contra¬ 
diction. Hence, also Case 2 cannot happen. 

Altogether, we have proved |i?i,[c, A]| = [S'!, [m, X], S'21, [c, A] |. Consider the 
path 

T := i?i, [c. A], A22, b,Z], [a;,Z],r(A'), [«,X]. 

By construction, |T| = k. Since level([M,X]) > k, the path T cannot be strongly 
simple. By construction, the paths Ri and ^(S'^) are not strongly disjoint. Let 
[d, A] be the first node on Ri such that one of the nodes [d, A] or [d, B] is on 
r(S'2). If [d, A] is on r(S'2) then a strongly simple path from s to [u, X] shorter than 
level([M,X]) can easily be constructed, a contradiction. Hence, [d,B] is a node on 
r(S'2) and hence, [d. A] is on S'2. Let Ri = Rn, [d, A],i?i2 and S'2 = S'21, [d, A],S'22. 
Consider the path 


T := Ru, [d. A], ^^2, [x, Z], [y, Z],r{S2), [u, X]. 

By construction, the path T is strongly simple. Since S' is a shortest strongly simple 
path from s to [x, Z] it holds that |i?ii, [d, A]| > |S'j, [u, X], S'21, ^]l- I’^ same 

way as we have proved |i?i, [c, A]| = [S'!, [m,X], S'21, [c, A]|, we prove |i?ii, [d, A]| = 
|S'j, [m, X], S'21, [d. A] I • Hence, |T| = k. This contradicts our assumption that 
level([-u,X]) > k. Altogether, the first property is proved. 

Assume that the second property is not fulhlled. Let m G H be a node with 
smallest second level such that the second property is not fulhlled with respect to 
the node u. This means that 

1 . I(level([M, A]) -|- level([M, B]) — 1 ) = I and there is a shortest strongly simple 
path P = P', [n,X], [m, X] with level2(M) = level([M,X]) but 

2 . P has not been constructed after the termination of Part 2 of Phase / 1 or 

level2(M) has not been computed correctly. 
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With respect to level([w,X]), the situations level([u,X]) = level([u,X]) — 1 and 
level([n,X]) < level([n,X]) — 1 are possible. We shall investigate both cases one 
after the other. 

Case 1: level([n,X]) = level([M,X]) — 1. 

Two subcases are to consider, level([n,X]) = leveli(n) and level([n,X]) = 
level2(n). 

Subcase 1.1: level([n,X]) = leveli(n). 

First we shall prove that level([n,X]) = level([M,X]) = 1. If level([n,X]) > I 
then level([M, X]) = level([n, X]) + 1 and level([n, A]) + level([M, B]) = 21 + 1 imply 
level([M, X]) < /. Let R be any shortest strongly simple path from s to [u, X]. Since 
leveli(n) > /, the nodes [n,A] and [n, 5] cannot be on R. Hence, i?, [n,X] would 
be a strongly simple path from s to [n,X] shorter than leveli(n), a contradiction. 
Fnrthermore, level 2 (n) = level([n,X]) + 1 and level([n,H]) + level([n, H]) = 2/ + 

1 imply level([n,X]) = I and also level([M,X]) = 1. Hence, dnring Part 1 of 
Phase I + 1, Case 2 applies such that the pair [n,X], [n,X] is inserted into E{1). 
Therefore, during Part 2 of Phase / + 1, the backward search with respect to the 
pair [u, X], [u,X] has been performed. 

By the induction hypothesis, at the beginning of Part 2 of Phase /+!, all shortest 
strongly simple paths from s to [u, X] and all shortest strongly simple paths from s 
to [u, X] are constrncted. Since dnring Part 2 of Phase / + 1, the backward search 
with respect to the pair [u, X], [u,X] is performed, the edges ([u,X], [u, X]) and 
([u,X], [u,X]) are inserted into the layered graph. Hence, after the termination 
of Part 2 of Phase I + 1, the path P = s,..., [u,X], [u,X] is constructed. By 
construction of the algorithm MBPS, level([u,X]) has been computed correctly as 
well. This contradicts onr assnmption such that Subcase 1.1 cannot happen. 

Subcase 1.2: level([u,X]) = level 2 (u). 

Let R = /?', [m,X] be any shortest strongly simple path from s to [m,X]. Con¬ 
sider the path Q := i?', [u, X], [u,X]. 

If Q is a shortest strongly simple path from s to [u,X] then level([u,H]) -|- 
level([u, H]) = 2/ -|- 1. By the choice of the node u, all shortest strongly simple 
paths from s to [u,X] has been constructed. By the induction hypothesis, R is 
also constructed. Hence, during the backward search which constructs the path 
s, P', [u,X], the edge ([u,X], [u,X]) is added to the layered network. Hence, the 
path P is also constrncted, a contradiction. Therefore, Q is not a shortest strongly 
simple path from s to [u,X]. 

If Q is strongly simple then level([u, X]) < |(5|. Hence, level([u, H])-|-level([u, P]) < 
level([M, H]) -|- level([u, B]) = 21 + 1. By the induction hypothesis, level([u, X]) has 
been compnted and all shortest strongly simple paths from s to [u,X] have been 
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constructed. Furthermore, leveli(M) = level([M,X]) < /. Hence, by the induction 
hypothesis, all shortest strongly simple paths from s to [u, X] have been con¬ 
structed and the edge ([m,X], [n,X]) has been considered during Part 1 of Phase 
level([n, X]) -|- 1. Therefore, during Part 2 of a Phase k where k < I + 1, the back¬ 
ward search has been performed with respect to the pair [n, X], [u, X] and the edge 
([n, X], [u, X]) has been inserted into the layered network such that the path P has 
been constructed after the termination of Part 2 of Phase k, a contradiction. 

Hence, the path Q = R', [m,X], [n,X] is not strongly simple. Since R', [n, X] is 
strongly simple, the node [n,X] or the node [n,X] has to be on the path R'. Let 


r R[,[v,X],R'^ if[v,X]eR' 
\ i?;,[n,X],i?' if[n,X]ei?'. 


If [n,X] is on R' then level([n,X]) < |(5|. Exactly as in the case “Q is strongly 
simple but IQI > level([n, X])”, we prove that the path P has been constructed 
after the termination of Part 2 of Phase I -|- 1 getting a contradiction. Hence, 
[n,X] is on R'. But then, R[, [n,X], [n, X] would be a strongly simple path from 
s to [n,X] shorter than level([M,X]), a contradiction. Therefore, Subcase 1.2 and 
hence. Case 1 cannot happen. 

Case 2: level([n,X]) < level([M,X]) — 1. 

Let Q, [uiX] be any shortest strongly simple path from s to [m,X]. Let [x,B] 
be the last node on P', [n,X] such that the length of the subpath from s to [x, B] 
of P is equal to level([a;, P]). Let P = Pi, [x, P], [y,A],P 2 , [n,X], [m,X], let R := 
[u,X], [v,X],riP 2 ), [y,B] and let S := Q, [u,X], [n,X],r(P 2 ), [y,B]. 

First we shall show that for all nodes [z.,Y] on P, the level of [z, Y] is not larger 
than the length of the subpath S'([ 2 :,y]) := Q,[u, X],r{P 2 i),[z,Y] of S. Then 
we shall use this fact to prove that all edges on the path R are contained in the 
current layered graph. Since level([|/, P]) < |S'([|/, P])|, the pair [x, B],[y, B] has 
been considered during Part 2 of a Phase k where k < I + 1. Hence, the path P is 
constructed after the execution of Part 2 of Phase k, a contradiction. 

To prove level([ 2 ;, H]) < |S'([ 2 ;, X])! assume that level([ 2 :, X]) > |S'([ 2 :,X])| for 
a node [z,Y] on R. Then S'([ 2 :,X]) cannot be strongly simple such that Q and 
P' ;= r(P 2 i), [z,Y] are not strongly disjoint. Let [p, H] be the hrst node on Q such 
that one of the nodes [p, A] or [p, P] is a node on P'. Let 


Q = Qi: [P: Q 2 and P' 


P{,[p,A],P' if[p,H]eP' 
Pj,[p,P],P' if[p,P]eP'. 


If [p, H] G P' then T := (5i,[p, H],P 2 would be a strongly simple path from s 
to [ZiY] shorter than |S'([ 2 :, X])|, a contradiction. Hence, [p, P] G P'. But then, 
T := Qi, [p, H], r(Pj), [n, X] would be a strongly simple path from s to [n, X] shorter 
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than level([M, X]), a contradiction. Altogether, we have proved that level([;s, F]) < 

|s(kr])|. 

Assume that an edge {[z, F], [z', F]) is not contained in the current layered net¬ 
work after the execution of Part 2 of Phase / + 1. Since level([z, F]) < |S'([z, F])|, 
we obtain level([F, F])+ level([ 2 ;, F]) < level([M, X])+ level([M, X]) — 1 = 21. Hence, 
^(level([F, F]) + level([ 2 :, F])) < 1. Therefore, the pair [F, Y][z, Y] has been consid¬ 
ered and the edge {[z^Y], [z',Y]) has been inserted into the current layered graph 
during Part 2 of a Phase k where fc < / + 1, a contradiction. Altogether we have 
shown that the path P is constructed after the execution of Part 2 of Phase / 1. 

Moreover, since the pair [x^B\, [y,B\ has been considered during Part 2 of Phase 
/ 1, level([M, X]) has been computed correctly because of the consideration of this 

pair of nodes. This shows that Case 2 cannot happen as well. This proves the 
second property. 

Altogether, the theorem is proved. | 

3.3 An implementation of an entire phase 

First we shall describe the implementation of MBFS, and then we shall show how 
to combine MBFS and MDFS to get an implementation of an entire phase. 

Obviously, Part 1 of all phases can be implemented in such a way that the total 
time is bounded by 0{m -|- n). For the implementation of Part 2 of all phases, we 
have to describe how to implement the search on the back-paths, starting in [n, Z\ 
and [w,F], until DOM({[n,F], [ta, Z]}) is reached. Most importantly, since we do 
not know DOM({[n, Z], [ta, Z]}) in advance, meaning that DOM({[n, Z], [w, Z]}) 
has to be computed simultaneously, we have to take care that the search does 
not continue beyond DOM({[n, Z], [w, Z]}). Note that by Lemmathe subpaths 
from s to DOM({[n, Z], [w, Z]}) are always shortest strongly simple paths from s 
to DOM{[n, Z], [w, Z]}. Hence, we can perform a breadth-hrst search on the back- 
paths until the current level contains exactly one node. By Lemma this node 
has to be DOM({[n, Z], [w, Z]}). 

With respect to the efficiency, at the moment when the search meets a node 
[n,X] for which level([M,X]) = level 2 (M) has been dehned, we have to compute 
efficiently the next node of the search having the property that its level is not 
dehned. By dehnition, this node is DOM([n,X]). As a consequence of Lemma 
we can maintain these nodes by a data structure for disjoint set union such that for 
the computation of DOM([n, X]) one find operation would suffice. In that case, an 
extensible edge would be stored. Using disjoint set union [3l], the total time can 
be bounded to be 0{{m + n)a{m,n)) where a is the inverse Ackermann function. 
Using incremental tree set union [12], we obtain a total time bound of 0{m + n) 
for the computation of the next node such that its level is not dehned. 
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The levels of the nodes [n, Z] and [w, Z] have not to be equal. If the level of 
the two nodes are different then we start the breadth-hrst search at the node with 
larger level. We always continue the search at the nodes with largest level until 
all nodes in the front of the search are on the same level. But we do not need to 
perform a precise breadth-first search. Hence, we can organize the search in the 
following way: 

1. In the front of the search, continue the search always in a node [m, W] such 
that there is another node \p,Y] in the front of the search with level([p, H]) < 
level ([m, X]). 

2. If the front of the search contains exactly one node then stop the search. 

Next, we shall combine MBFS and MDFS for the implementation of an entire phase. 
Knowing Gm, a maximal set of up to s and t pairwise disjoint shortest strongly 
simple paths from s to f in Gm can be computed using MDFS in 0{m + n) time. 
Every time, a strongly simple path P from s to t is found, all nodes [n,H], [v,B] 
with [n, H] G P or [n, P] G P and all incident edges are deleted from Gm- If a node 
gets zero indegree or zero outdegree, then also this node, and all incident edges, 
are deleted. Altogether, we have obtained the following theorem. 

Theorem 6 A maximum matching in a general graph G = (V, E) can be computed 
in 0{^/n{m + n)) time and 0{m + n) space, where \V\ = n and \E\ = m. 

4 The primal-dual method for the 

maximum weighted matching problem 

Let G = (y,E) be an undirected graph. If we associate with each edge (i,j) G E 
a weight w{i,j) > 0 then we obtain a weighted undirected graph G = (y,E,w). 
The weight w{M) of a matching M is the sum of the weights of the edges in M. 
A matching M P E has maximum weight iff J2(i j)&M' G Xlp j)eM j) 

for all matchings M' C E. Given a weighted undirected graph G = {V,E,w), the 
maximum weighted matching problem is Ending a matching M C P of maximum 
weight. Our goal is to develop a method for the computation of a maximum 
weighted matching in a given weighted undirected graph. 

4.1 The description of the method 

Let G = (y,E,w) be a weighted undirected graph. Let E = {Ei, E 2 ,..., E^}, 
Pj C P be a family of pairwise distinct subsets of P. With each node i G K 
we associate a node weight 7r(i) > 0. Furthermore, with each edge set Ei G E, we 
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associate a set weight iJ,{Ei) > 0. These new variables are called dual variables. The 
values of the dual variables are treated in such a way that the following invariant 
is always maintained. 

• w{i,j) < 7r{i) + 7r(j) + E(*j)eE, ^r all {i,j) G E. 

The right side of this inequality is the dual weight d{i,j) of the edge {i,j) E E. We 
dehne the dual weight d{M) of a matching M by 

d{M) := d{i,j). 

Note that because of the invariant, always w{M) < d{M) for all matchings M E. 
With respect to an arbitrary matching MCE, the maximal possible contribution 
of the node weight n{i) to d{M) is n{i) since i is adjacent to at most one edge in 
M. Note that \Ei fl M\ < c{Ei) where c{Ei) is the size of a maximum cardinality 
matching with respect to Ei. Hence, the maximal possible contribution of the set 
weight fJ,{Ei) to d{M) is c{Ei)fj,{Ei). Hence, 

^7r(i) + Y c{Ei)fi{Ei) 

i£V Ei&T 

is always an upper bound for the dual weight of any matching of G. Therefore, 
with respect to a matching M, 

w{M) = ^7r(f) + Y c{Ei)^i{Ei) 

i&V EieT 

implies that the matching M has maximum weight. The question is now: When 
with respect to a matching M, this equation holds? 

The right side of the equation contains for each edge in M its dual weight. Since 
the dual weight of an edge is at least as large as its weight, we obtain the necessary 
condition d{i,j) = w{i,j) for all edges {i,i) G M. Since all summands in both sums 
of the right side of the equation are non-negative, the node weight n{i) has to be 
zero for all M-free nodes i E V. Furthermore, for all Ei such that \Mr\Ei\ < c{Ei), 
the set weight n{Ei) has also to be zero. On the other hand, these conditions are 
fulhlled with respect to a matching M for which w{M) = d{M). Altogether, we 
have obtained the following necessary and sufficient optimality conditions: 

1. r{i,j) := d{i,j) - w{i,i) = 0 for all (i, j) G M, 

2. Ti{i) = 0 for all M-free nodes i E V, and 

3. fv{Ei) > 0 for El E E implies \Ei 0 M\ = c{Ei). 


38 



The value r(i, j) is called the reduced cost of the edge {i, j). Because of the invariant 
which is maintained, the reduced cost of an edge is always non-negative. 

The primal-dual method for the weighted matching problem can be separated 
into rounds. The input of every round will be a matching M, a set of pairwise 
distinct subsets of E such that > 0 for all Ei G E, and values for the dual 

variables which fulhll the hrst and third optimality conditions with respect to the 
matching M. The second optimality condition can be unsatisfied with respect to 
some free nodes. Our goal within a round is to modify M and the values of the 
dual variables such that Conditions 1 and 3 remain valid and the number of nodes 
violating Condition 2 is strictly decreased. 

A round divides into two steps, the search step and the extension step. The 
search step tries to improve the current matching by hnding an augmenting path 
P such that the number of free nodes with node weight larger than zero can be 
decreased by the augmentation of P. Since Condition 1 has to be maintained, 
the search step can only be performed on edges with reduced cost zero. If the 
search step cannot improve the current matching then the extension step changes 
the values of some dual variables using an appropriate value S. The extension step 
can decrease the reduced cost of some edges to zero. Hence, it is possible that the 
next search step will hnd an augmenting path. 

During the search step, we use MDFS. Hence, we dehne with respect to the 
current matching M the weighted directed bipartite graph Gm '■= (y',EM,w) 
where 

V' := {[v, A],[v, B] \ V E V} U {s,t} s,t^V,Sy^t 
Em ■= {i[v,A],[w,B]),{[w,A],[v,B])\{v,w)EM} 

^ {{[x, B],[y, A]), {[y, B],[x, A]) \ {x,y) eE\M} 

U {(s, [u, B]), {[v, A],t) \ V E V is M-free} . 

This means that Em contains for each edge in E two copies. Both copies of the edge 
{i,j) obtain weight w{i,j), dual weight d{i,j), and hence, reduced cost r{i,j)- We 
arrange that edges with tail s or head t have always reduced cost zero. According 
to Condition 1, it is only allowed to consider augmenting paths where all edges on 
these paths have reduced cost zero. Hence, the input graph G*j^ = (y',Elj,w) of 
the search step will be the subgraph of Gm containing exactly those edges in Em 
having reduced cost zero; i.e., 

El, := {([^,X], [j,X]) E Em \ r(*, j) = 0} U {(s, [t,B]), {[t,A],t) \ t M-free}. 

By dehnition, for each edge {i,j) E M it holds that {[i,A], [j,B]), {[j,A], [i,B]) E 
El,. We start with the empty matching 0 and define the graph G(d = (y',E(D,w) 
as described above. Let W := We initialize all node weights 7i{i) 
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by Altogether, we obtain the input graph G*^ = {V',E^,w) for the first search 
step where 

E;= {{[t,B],[j,A]),{[j,B],[t,A])\{[t,B],[j,A])eE^, w{t,j) = W} 

U{(s, {[i,A],t) \i^V}. 

At the beginning, the family E of subsets of E will be empty such that no set weight 
has to be defined. During the execution of the algorithm, the needed elements of 
E and the corresponding set weights will be defined. Whenever to some edge set 
El ^ E a value fi{Ei) > 0 is associated, Ei enters E. As soon as iJ,{Ei) gets to be 
zero, the set Ei is removed from E. 

A search step terminates with a matching M, a weighted directed graph Gm = 
(!/', Emi w), a current subgraph G*^ = (W, E\j, w) of Gm such that no M-augmenting 
path P is contained in G\j, a current set E where the set weight of each set in E is 
positive, and values for the dual variables. This is the input of the next extension 
step. 

For the treatment of the extension step consider the expanded MDFS-tree Te^p, 
computed by the last MDFS on G\j. Note that this MDFS was unsuccessful; i.e., 
no path from the start node s to the target node t has been found. The goal is to 
add edges to Texp such that possibly an augmenting path can be found. Therefore, 
we have to decrease the reduced cost of edges with positive reduced cost. Such an 
edge (hi) has to be in i? \ M and [i, B] has to be in Texp- Let 

Va-.= {[i,A]\ieV} and Vb ■-= {[i,B]\i eV}. 

Furthermore, let 

At := Va a Texp and for a moment Bt := Fb H Texp- 

Later, we shall see that according to the optimality conditions which we have to 
maintain, some nodes in Fg fl Texp are not allowed to be a node in Bt- Hence, we 
shall modify the definition of Bt by removing exactly these nodes. 

The idea is to decrease the reduced cost r{i,j) of all edges {i,j) with positive 
reduced cost and [i, B] G Bt by an appropriate value S. This is done by decreasing 
the node weight n{i) by S for all nodes i with [i, B] G Bt- As a consequence of the 
decrease of the node weight the reduced cost of each edge e in G*m with end 
node [i, A] or [i, B] becomes negative. Because of the invariant maintained by the 
method, such a reduced cost has to be increased until its value is zero again. We 
distinguish two cases: 

1. The other end node of e is in At but not in Bt- 

2. The other end node of e is in Bt- 
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If Case 1 is fulfilled then we can increase the reduced cost of the edge e by increasing 
the node weight of the other end node of e by 5; i.e., we increase 7r{j) by 6 for all 
nodes j such that [j, A] G At and [j, B] ^ Bt- Note that increasing the node weight 
7r(j) of a node j increases the reduced cost of each edge (f, j). If (i, j) E E \ M 
then r{i,j) has only to be non-negative. Since r{i,j) was non-negative before the 
increase of vr(j), it is also non-negative after the increase. If {i,j) G M then r{i,j) 
has to be zero. Note that [j, A] G At implies for the unique node [f , B] with 
(f, j) G M that [i, B] G Bt- Hence, 7i{i) has been decreased by the same value such 
that r{i,j) = 0 before the change of the dual variables implies that r{i,j) = 0 after 
the change. 

If Case 2 is fulhlled then the reduced cost r{i,j) is decreased by 2S. This can 
be corrected by increasing the set weight fi{Ei) of exactly one set Ei containing the 
edge (i, j) by 25. Two questions have to be answered. 

1. What is the accurate edge set Ei for increasing its set weight? 

2. What is the appropriate value 52 

To answer the hrst question, let us consider MDFS which is used as a subroutine 
during the search step. Review the dehnitions and properties of the sets and 

D[q,A] as given on Pages and respectively. At the moment when an edge set 
Eg is chosen to obtain a positive set weight fi{Eq), the set Eg corresponds to the 
current set T*[g,A] as defined on Page [T7| and Eg enters E. Note that E can contain 
edge sets which are generated during a previous run of MDFS. The corresponding 
set has not to be current with respect to the current run of MDFS. Hence, 

we introduce an analogous terminology to “current” with respect to the edge sets 
in E. A set Eg E E is called maximal iff Eg ^ Eg^ for all Eg/ E E \ {Eg}. Now, it 
is useful to investigate the structure of a set Dig^A]- Let 

Diq,A] ■= {q}A{pEV\ \p, A] G D[q^A]} 


and 

D[qA] [p^ B]\pE D[g^^^}. 

The node q is the unique node p E D'^qA] such that p is end node of an edge 
(r,p) G M with r ^ B'^qA]- ^ path P runs through a set then there is an 

edge such that P enters D^^g^A] using this edge and also an edge such that P leaves 
B[q,A] using that edge. We say that a path P enters or leaves D^g^A] via an edge in 
E\M ii the used edge ([x, R], [p. A]) corresponds to an edge (x, y) E E \ M. Let 
(r, g) G M be the unique matched edge with end node q. During the execution of 
MDFS, for an M-augmenting path P, there are three possibilities to run through 
a set D[q,A]- 
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1. P enters via the matched edge ([r, A], [q,B]) and leaves -D[q,A] via an 

edge in £■ \ M. 

2. P enters via an edge in E\M and leaves -D[gyi] via the matched edge 

{[q,A],[r, B]). 

3. P enters and leaves Dig^A] via an edge in \ M. 


If an M-alternating path R enters P>[q,A] via the edge ([r, ^4], [q, B]) then, by Lemma 
3, for all n G -Djg y^] , [v,B] E Bt- Therefore, with respect to each edge ([i, X], [j, X]) 

in E\j with both end nodes in -D[g,A], the node weights 7r(z) and 7r(j) have been 
decreased by 5 snch that we have to increase the set weight of exactly one edge set 
containing the edge (i,j) by 25. Note that for all n G Id there exists at most one 
cnrrent D\^q^A] snch that v G Hence, we dehne the edge set Eq corresponding 

to D[^,a] where -D[g,yi] cnrrent by 


Eq X n E. 

If we have to increase the set weight with respect to an edge ([h-B], ^ E\j 

then we choose the edge set Eq where -D[q,A] is the current set such that i,j G -Dj^ y^^. 
Hence, for all current sets Dyq^A], we perform the following operation: 

• If Eq is not already in E then Eq is inserted into E with fJ^{Eq) := 26. 
Otherwise, fi{Eg) is increased by 26. 

The following lemma shows that in this situation, the current set -D[q,yi] always 
exists. 


Lemma 7 Let [i,B],[j,B] G Hr and {[i, B],[j, A]), ([j, B],[i, A]) G Then 

there exists a current set D^^g^A] such that i,j G 

Proof: Note that {[i, B],[j, A\),{[j, B],[i, A]) G El^ iff r{i,j) = 0. With respect 
to the positions of the nodes [i, B] and [j, B] in the MDFS-tree T, exactly the 
following two situations are possible: 

1. There is no path P from the root s to a leaf of T such that both nodes are 
on P. 


2. There is a path P from the root s to a leaf of T such that both nodes are on 
P. 


If the hrst situation arises then there exists a unique node in T such that both 
nodes [i, B] and [j, B] are successors of this node in T and no successor of this node 
in T has this property. Since no M-augmenting path has been found during the 
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last MDFS, this node cannot be the node s. Hence, this node has to be a node 
[u,B] G Vb- Let s,Q, [u,B] be the path from s to the node [u,B] in T. Further¬ 
more, let [u, -B], -Pi, [i, B] and [u, B], P 2 , [j, B] be the paths in T from [u, B] to [i, B] 
and [j,-B], respectively. By the dehnition of and because of the paths 

s,Q, [m,-B],Pi, [i,-B], [j,H],r(P 2 ), [u,A] ands,Q, [m,P],P 2 , [j,P], [i,H],r(Pi), 

= ['n,H] if [m, H] ^ At. If [u,A\ G Ht, i.e., PUSH([m,H]) has been 
performed, then, by Lemma ~ ~ L]^u,a]- By the dehnitions of the 

current set P[g,A] and the set P|^ , the assertion follows directly. 

Assume that the second situation arises. W.l.o.g., let [j, P] be a successor of 
[i, B] in T. Let s, Q, [i, P], P, [j, P] be the path from s to [j, P] in T. 

If [i, A] ^ Ar then, because of the path s, Q, [i, P], [j. A], r(P), [i. A] and the 
dehnition of -Lp^A], it follows -L[j;a] = [*,A]. By the dehnition of P[i,A], the node 
[j. A] is contained P[i,A]- The dehnition of P{j^] implies i,j G Pj^A] such that the 
assertion follows directly. 

If [i,A] G At then, by Lemma L^^a] = -L[i,A]- Hence, the assertion follows 
directly. | 


Let us examine the ehect of the augmentation of the path P to the number of 
edges in the current matching with both end nodes in P|^ . If P enters and leaves 

P[g,A] via an edge in P \ M, then this number decreases by one. In the other two 
cases, this number does not change. Hence, the augmentation of an augmenting 
path which goes through D^q using the unique matched edge with one end node in 
and the other end node not in P|^ is always allowed. But the augmentation 

of an augmenting path which enters and leaves -D[g,A] via an edge in P \ M is only 
allowed if fi{Eq) = 0; i.e.. Eg is not contained in P. Now, we shall determine the 
accurate value for 6. 

Since all node weights have to be non-negative, 6 cannot be larger than the 
node weight of any M-free node i. Note that with respect to an M-free node i, 
always [i,B] G Bt is fulhlled. Since all nodes are initialized with the same node 
weight and always the node weights of all nodes in Bt are decreased using the same 
value, all free nodes have the same node weight and no non-free node has a smaller 
node weight. Hence, 6 < 7r(z) where i is M-free implies that after the change of 
the dual values, all node weights are non-negative. 

Let [i, P] G Bt- Then 7i{i) will be decreased by 6. For {i,j) G P, in dependence 
of the status of the other end node j with respect to Bt and A^, three situations 
can arise: 

If [j,B] ^ Bt and [j. A] ^ At then n{j) will not be changed by the extension 
step. Hence, we have to choose 6 := r{i,j) to decrease the reduced cost r{i,j) of 
the edge (i, j) to zero. 

If [j,P] ^ Bt and [j. A] G At then 7r(j) will be increased by S. Hence, inde¬ 
pendently of the choice of 6, the reduced cost r{i,j) does not change. 
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If [j, B] G Bt then 7r(j) will be also decreased by S. Hence, we have to choose 
S := |r(h j) to decrease r{i,j) to zero. 

Note that 6 has to be chosen in such a way that after the extension step r{i, j) > 
0 for all edges {i,j) G E. Hence, 5 should not be larger than the minimal reduced 
cost with respect to edges (i,j) with [i,B] G Bt, [j,B] ^ Bt and [j, H] ^ At, and 
also not larger than the half of the minimal reduced cost with respect to edges 
{i,j) with [i,B], [j,B] G Bt, r{i,j) > 0, and there is no current a] such that 
Note that i,j G with respect to a current set implies 

that r{i,j) does not change since jJi{Eq) is increased by 25. Because of Lemma 
r{i,j) = 0 implies that there is a current £*[ 5 ,^ 4 ] such that i,j G D'^qA]- 

Since MDFS obtains in addition to the input graph the current set E and 
the corresponding set weights, because of the maintenance of the third optimality 
condition, the following holds true with respect to a maximal Eq & E\ 

Assume that during the MDFS no path uses the edge {\r,A\, [q,B]) to enter 
D\^q^A]', be., [q,B] ^ Tf.xp- If there is a path R entering D[g,A] via an edge in E\M 
then fJ,{Eq) > 0 implies that R has to leave E>iq,A] via the edge {[q, A],[r, B]), 
independently if the node [q,A] is already pushed or not. Since (g,r) G M and 
therefore {[q,A], [r,B]) G El^, this is always possible. Note that more than one 
such a path can run through D^q^A]- Depending on the nodes in D^^g^A] which are 
entering nodes of such a path R, the following situations can happen with respect 
to a node v G D'^gA]- 

a) [v,B] G Bt, 

b) [v,B] ^ Bt but [v,A\ G At, or 

c) [v, B] ^ Bt and [v. A] ^ At- 

The problem to solve is the following: How to change the node weights of the nodes 
in Each entering node [v,A\ of D]^q^A] is the head of an edge {[x,B], [n, A]) 

where [x, B] G Bt- Hence, n^x) is decreased by 5. According to the hrst optimality 
condition, 7r(n) has to be increased by 5. Possibly, there are edges in Eg with exactly 
one end node is an entering node, with both end nodes are entering nodes or with 
no end node is an entering node. With respect to all these cases, the node weights 
and fi{Eg) have to be changed in such a way that the optimality conditions remain 
to be valid. For doing this, we increase 7i{v) by 6 for all v G D'^gA]- Since we have 
increased the reduced cost of all edges in Eg by 25, we decrease fi{Eq) by 25. Since 
n{Eq) has to be non-negative, 5 has to be chosen such that before the change of 
the dual variables, n{Eq) > 25. 

Remember that we consider the situation that [q, B\ ^ T^xp- Possibly, there 
exist nodes [i, B] G D\^g^A] which are in T^xp- Since we cannot decrease 7r(i) by 5 
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and increase 7 i{i) by 6 at the same time, it is not allowed for such a node to be in 
Bt- Hence, we dehne 

Bf := {[i,B] I [i,B] e D[g^A], Eg maximal and [q,B] ^ Tf,xp} 
and we redehne 

Bt := {Vb H T^xp) \ Bf. 

Altogether, 6 can be dehned in the following way. 


(5o := '7r(i), where i is M-free, 

(5i := min{r(b j) | [j, A] ^ At, [j, B] ^ Bt and [i, B] e Bt], 

62 := min{r(bj) | [i,B], [j,B] G BT,r{i,j) > 0 and there is no current 

L>[q,A] with i,i e D[q jf^], and 

^3 := mm{fi{Eq) \ Eg G E maximal, [q,B] ^ Bt and [q,A] G At}. 


Then we dehne 


6 := min 





Altogether, we have obtained the following extension step. 


(5o := n{i) for an M-free node i; 

:= min{r(i, j) | [j,A] ^ At, [j,B] ^ Bt and [i,B] G Bt}', 

82 := min{r(i,j) | [i,B], [i,B] G BT,r{i,j) > 0 and there is no current 
Dyg^A] with i, j G 

53 := min{/i(i?q) \ Eg ^ E maximal, [q,B] ^ Bt and [q,A 8 \ G At}', 

5 := min{5o,(5i, f, f }; 
for all [i, B] G Bt 

do 

■n{i) := Ti{i) — 5 

od; 

for all {[i,B] ^ Bt, [i. A] G At) or {i G D'^gAY Eg E E maximal with 

[q, B] ^ Bt but [q. A] G At) 

do 

Ti{i) := Ti{i) -|- 5 

od; 

for all current T*[5 ,a] and [q, B] G Bt 

do 

•= + 25 

od; 

for all maximal Eg G E, [q, B] ^ Bt and [q. A] G At 
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do 


fi{Eg) := ^{Eg) - 26 

od. 

The correctness of the described primal-dual method follows from the discussion 
done during the development of the method. Note that the primal-dual method 
described above is equivalent to the method developed by Edmonds [T3] . 

4.2 An implementation of the primal-dual method 

First we shall determine an upper bound for the number of dual changes which 
can occur between two augmentations in the worst case. Then we shall describe 
the implementation of the search steps between two augmentations. Finally, we 
shall describe an implementation of the computation of 6 and the update of the 
dual variables. With respect to the determination of the upper bound, we have to 
consider four cases. 

Case 1: 6 = 6 o 

Then after the change of the dual variables, n{i) = 0 for all M-free nodes 
i G F. Therefore, the current matching M is of maximum weight and the algorithm 
terminates. Hence, Case 1 occurs at most once. 

Case 2: 5 = (5i 

Then during the next search step, at least one new node [j, A] enters Ax. Hence, 
Case 2 occurs at most n times. 

Case 3: 6 = ^ 

Then during the next search step, at least one new edge (z, j) enters E*j^. Fur¬ 
thermore, at the moment of the definition of the value 5, there is no current Dyg^A] 
with i,j G D'^gA]- Lemma ItI shows that after the next search step there exists a 
current set such that^j^^j contains both nodes i and j. This means, by the 

union of two smaller current sets, a larger current set is obtained. Between two 
augmentations, the number of such unions is bounded by n — 1. Hence, Case 3 
occurs at most n — 1 times. 

Case 4-' 6 = Y 

Then at least one edge set Eg leaves the family E. As long as [q, B] ^ Bt, no 
edge set Eg can enter the family E. But if [q, B] is pushed, a corresponding edge 
set Eg cannot contribute to the definition of the value 63 before [q, B] leaves Bt 
again. This cannot happen before the next augmentation. Hence, Case 4 occurs 
at most n times. 
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Altogether we have shown that the number of dual changes between two aug¬ 
mentations is bounded by 3n. 

With respect to the implementation of the search steps between two augmenta¬ 
tions, MDFS has to be adapted to the situation that the current search path enters 
the set -D[g,yi] via an edge m. E\M with respect to a maximal edge set Eq G E. 
According to the discussion during the development of the primal-dual method, we 
have to jump to the node [q, A] and to continue the search using the unique edge 
([g. A], [r, B]) G E\j. This can be organized using a data structure for disjoint set 
union. But this data structure has also to support a further operation. 

If according to an extension step, n{Eq) becomes zero, the corresponding edge 
set leaves the family E. Note that Eq is a maximal set in E. If there is a set 
Eq/ G Eq in E then another set in E becomes to be maximal. This means that we 
have to undo the union operation performed with respect to the set Eq. 

With respect to the implementation of MDFS, we have introduced a data struc¬ 
ture for disjoint set union which uses the weighted union heuristic. Our goal is now 
to extend this data structure to support also deunion operations such that the time 
used for the deunion operations will be, up to a small constant factor, the same as 
the time used for the union operations and each find operation uses only constant 
time. This can be done in the following way: 

During the execution of an union operation, instead of changing a pointer, we 
add a new pointer. The current pointer of an element will be always the last 
created pointer. Since the pointer of any element is changed at most logn times, 
for each element, at most logn extra pointer are used. The time used for the 
union operations remains essentially the same. A deunion can be performed by the 
deletion of the current pointers created during the corresponding union operation, 
the update of the set sizes and the update of the name of the larger subset. The time 
used for a find operation remains to be constant. Altogether, the data structure 
for union-find-deunion can be implemented such that the time used for at most m 
find, n union and n deunion operations is 0 {m + nlogn). 

After an extension step, the last MDFS can be continued instead of to start a 
new MDFS. Next we shall describe an implementation of the computation of the 
value S and the update of the dual variables. 

Since all M-free nodes i have the same dual weight, (5o can be computed by the 
consideration of any M-free node in constant time. 

For the computation of (5i, we maintain a priority queue Vi which contains for 
all [j, A] ^ At with [j, B] ^ Bt an item which points to a list containing exactly 
the edges {i,j) with [i,B] G Bt of minimum reduced cost within all such edges, 
if such an edge exists. The associated key with this item is the minimum reduced 
cost of these edges. Furthermore, using an array of size n, we have direct access 
to the list in Vi containing all edges with end node [j. A], if such an edge in Vi 
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exists. Note that each extension step decreases the keys of all items in the priority 
qneue by the cnrrent 6 . It is usefnl to maintain the property that always the key 
of all elements in Vi has to be decreased by the same amonnt. Hence, we maintain 
the sum Ai of all dual changes done so far. If we add an edge {i,j) to the priority 
queue then we dehne the modihed reduced cost of the edge to be r(i,j) + Ai. At 
the moment when a node [i, B] enters Bt, it is possible that we have to insert some 
edges into Vi. Hence, we update Vi with respect to [i, B] G B^ at that moment in 
the following way: 

For all edges {i,j) with [j,B] ^ Bt and [j, A] ^ At perform the following 
update operations: 

(1) If no element with respect to j is contained in the priority queue then we 
create a list which contains the element {i,j) with modihed reduced cost 
r{i,j) + Ai. We insert a new item which points to this list with associated 
key r{i,j) + Ai into Vi. 

(2) If Vi contains an item which points to a list containing edges with end node 
[j. A] with larger associated key than r{i,j) + Ai then delete all edges from 
the list, insert the edge by ([i, B], [j, A])) into the list and decrease the key of 
the item such that its value becomes r{i,j) + Ai. 

(3) If Vi contains an item which points to a list containing edges with end node 
[j, A] with associated key r{i,j) + Ai then add the edge ([f, B], [j, A])) to the 

list. 

(4) In all other cases do nothing. 

(5) If Vi contains some edges with respect to the node [i. A] then delete the 
corresponding item and the list of edges to which this item points. 

The number of deletions performed in Step 5 is bounded by the number of nodes in 
Va- If 5 = (fi, we have to delete at least one minimum key from Vi. Altogether, the 
total number of deletions is bounded by n. For the implementation of the priority 
queue Vi we can use Fibonacci heaps [H] or strict Fibonacci heaps ra such that 
with respect to the computation of all 5i’s between two augmentations, the used 
time is 0 {m + n log n). 

For the computation of all S 2 A, we maintain a priority queue V 2 which contains 
all edges (i,j) such that [i, B],[j, B] G Bt and r{i,j) > 0. Again, we can use 
Fibonacci heaps or strict Fibonacci heaps for the realization of the priority queue. 
Note that each extension step decreases all weights of the elements {i,j) in V 2 with 
the property that there is no current such that i,j G D'^gA] where 6 is 

the value chosen for the current extension step. Analogously to the manipulation 
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of Vi, we maintain with respect to 1^2 the sum A 2 of all dual changes done so far 
with respect to edges in V 2 and modify the weights in the appropriate manner. We 
update V 2 with respect to [i, B] G Bt at the moment when [i, B] is added to Bt 
in the following way: 

• For all [j,B] G Bt with r{i,j) > 0 add the edge {i,j) with modihed reduced 
cost r{i,j) + A 2 to the priority queue V 2 - 

Since at most m edges are inserted, the used time is 0{m). 

For the computation of 52, we have to hnd the smallest element (f,j) in V 2 
which has the property that there exists no current such that G D'^qA]- 

We use the priority queue V 2 and perform the following procedure: 

(1) hndmin; 

(* Let (f, j) be the output of hndmin.*) 

(2) if there exists F*[q,A] current with G 

then 

deletemin; 

goto (1) 

else 

^2 := - 82 

fi. 

If 5 = ^2 then we have to delete at least one minimal element from the priority 
queue V 2 - The number of deletemin operations performed during all computations 
of 82 between two augmentations is bounded by the number of edges in E. Alto¬ 
gether, the number of deletions performed during the computation of 82 is 0 {m). 
Each deletion can be performed in O(logn) time. Hence, the total time used for 
deletions is O(mlogn). Altogether, with respect to the computation of all 52 ’s 
between two augmentations, the used time is 0(m logn). 

For the computation of all 83 s, we maintain a priority queue V 3 which contains 
for all maximal Eq E E with [q,B] ^ Bt and [g. A] G At the value fi{Eq). We 
use a heap, a Fibonacci heap or a strict Fibonacci heap for the realization of the 
priority queue V 3 . Each extension step decreases all keys of the elements in P 3 . 
The amount is two times the current 8 . Hence, we can use the value A 2 dehned 
above and modify the keys in the appropriate manner. We update V 3 before the 
computation of 8 . We have to insert for all [q, A] G Va such that Eq is maximal and 
[q, A] was inserted into Texp after the last dual change and the last augmentation, 
respectively but [q, B] ^ T^xp the value ^{Eq) + A 2 . We have to delete for all [q, B] 
which are inserted to T^xp after the last dual change and the last augmentation, 
respectively the corresponding value if in V 3 such a value exists. Since the number 
of dual changes between two augmentations is bounded by 3n and [q, B] G Texp 
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can only leave Texp becanse of an angmentation, the nnmber of insertions and 
deletions is at most 2n. If 5 = 63 , we have to delete at least one minimal element 
from Vs- As observed above, the nnmber of snch deletions is bonnded by n. If 
there exists some sets Eq/ C Eg in E, some other edge sets in E become to be 
maximal. Hence, we have to perform the dennion operation which corresponds 
to the nnion operation performed with respect to the constrnction of Eq and to 
insert the resnlting maximal sets Eqi of E with key fJ-{Eqi) + A 2 . As observed 
above, the total nnmber of snch insertions is also bonnded by n. Each deletion and 
each insertion can be performed in O(logn) time. Altogether, with respect to the 
compntation of all (Js’s between two angmentations, the nsed time is 0(nlog?7,). 

We have proved the following theorem. 

Theorem 7 The primal-dual method can be implemented such that its time com¬ 
plexity is 0{nm\ogn). 

To get an implementation of time complexity 0{nm\ogn), we can nse ordinary 
heaps instead of Fibonacci heaps or strict Fibonacci heaps. Using Fibonacci heaps 
or strict Fibonacci heaps, the valne hi can be compnted snch that the nsed time 
between two angmentations is 0(m + nlogn). To get an implementation of time 
complexity 0 {n{m + nlogn)) the only critical point is the compntation of the 
valnes 62 snch that the time nsed between two angmentations is bonnded by 
0{m + nlogn). Gabow [H] has presented an implementation of Edmonds maxi- 
mnm weighted matching algorithm which nses complicated data strnctnres at the 
hrst SODA. The stated time complexity is 0{n{m + nlogn)). 

The description of some recent programs which implements Edmonds’ maxi- 
mnm weighted matching algorithm can be fonnd in P, 1301 1^ . 

Acknowledgments: I thank Ross McConnell for pointing ont a mistake in the 
description of the algorithm in [S] and Ari Frennd for indicating some points in [7] 
which have to be clarified. These hints have cansed myself to revise my whole work 
on graph matching algorithms. I thank Mathias Hanptmann for a carefnl reading 
of a preliminary version of the paper. 


References 

[ 1 ] Aho A. V., Hopcroft J. E, and Ullman J. D., The Design and Analysis of 
Computer Algorithms, Addison-Wesley (1974), 187-189. 

[2] Balinski M. L., Labelling to obtain a maximnm matching, in Combinatorial 
Mathematics and its Applications (R. C. Bose and T. A. Dowling, eds.). Uni¬ 
versity of North Carolina Press, Chapel Hill (1969), 585-602. 


50 


[3] Bartnik G.W., Algorithmes de couplages dans les graphes, These Doctorat 3® 
cycle, Universite Paris VI (1978). 

[4] Berge C., Two theorems in graph theory, Proc. Nat. Acad. Sci. U.S.A., 43 
(1957), 842-844. 

[5] Blnm N., A new approach to maximnm matching in general graphs, 17th 
ICALP (1990), LNCS 443, 586-597. 

[6] Blnm N., A new approach to maximnm matching in general graphs. Research 
report, Universitat Bonn (1993). 

[7] Blnm N., Maximnm matching in general graphs withont explicit consideration 
of blossoms. Research report, Universitat Bonn (1999). 

[8] Blnm N., A simplihed realization of the Hopcroft-Karp approach to maximnm 
matching in general graphs. Research report, Universitat Bonn (1999). 

[9] Cook W., and Rohe A., Compnting minimnm-weight perfect matchings, IN¬ 
FORMS Journal on Computing 11 (1999), 138-148. 

[10] Brodal G. S., Lagogiannis G., and Tarian R. E., Strict Fibonacci heaps, AAth 
STOC (2012), 1177-1184. 

[11] Even S, and Kariv O., An 0(n^'®) algorithm for maximnm matching in general 
graphs, 16th FOCS (1975), 100-112. 

[12] Edmonds J., Paths, trees, and flowers, Canad. J. Math, 17 (1965), 449-467. 

[13] Edmonds J., Maximnm matching and a polyhedron with 0,1-vertices, J. Res. 
Nat. Bur. Standards 69 B (1965), 125-130. 

[14] Fredman M. L., and Tarjan R. E., Fibonacci heaps and their uses in improved 
network optimization algorithms, J. ACM 34 (1987),596-615. 

[15] Freund A., private communication 2013. 

[16] Gabow H. N., Implementations of algorithms for maximum matching on non- 
bipartite graphs. Doctoral thesis. Comp. Sci. Dept., Stanford Univ., Stanford, 
Calif. (1973). 

[17] Gabow H. N., An efficient implementation of Edmonds algorithm for maximum 
matching on Graph, J. ACM, 23 (1976), 221-234. 

[18] Gabow H. N., Data structures for weighted matching and nearest common 
ancestors with linking, 1st SODA (1990), 434-443. 


51 



[19] Gabow H. N., and Tarjan R. E., A linear-time algorithm for a special case of 
disjoint set nnion, J. Comput. Syst. Sci., 30 (1985), 209-221. 

[20] Gabow H. N., and Tarjan R. E., Faster scaling algorithms for general graph¬ 
matching problems, J. ACM 38 (1991), 815-853. 

[21] Galil Z., Efficient algorithms for finding maximum matching in graphs, Gom- 
puting Surveys, 1986, 23-38. 

[22] Galil Z., Micali S., and Gabow H. N., An 0{EV log V) algorithm for finding a 
maximal weighted matching in general graphs, SIAM J. Comput. 15 (1986), 
120-130. 

[23] Gondran M., and Minoux M., Craphs and Algorithms, Wiley & Sons, (1984). 

[24] Hopcroft J. E., and Karp R. M., An algorithm for maximum matching in 
bipartite graphs, SIAM J. Comput., 2 (1973), 225-231. 

[25] Kariv O., An algorithm for maximum matching in general graphs, 

Ph.D. thesis. Dept, of Applied Mathematics, Weizmann Institute of Science, 
Rehovort, Israel (1976). 

[26] Kolmogorov V., Blossom V: A new implementation of a minimum cost perfect 
matching algorithm. Mathematical Programming Computation 1 (2009), 43- 
67. 

[27] Lawler E.: Combinatorial Optimization, Networks and Matroids, Holt, Rine¬ 
hart and Winston, (1976). 

[28] Lovasz L., and Plummer M. D., Matching Theory, North-Holland Mathematics 
Studies 121, North-Holland, New York (1986). 

[29] McGonnell R., private communication 2011. 

[30] Mehlhorn K., and Schafer G., Implementation of 0{nm logn) weighted match¬ 
ings in general graphs: The power of data structures, ACM Journal of Exper¬ 
imental Algorithms 7 (2002). 

[31] Micali S., and Vazirani V. V., An 0(a/|K| ■ |E|) algorithm for finding maximum 
matching in general graphs, 21st FOCS (1980), 12-27. 

[32] Peterson P. A., and Loui M. G., The general matching algorithm of Micali and 
Vazirani, Algorithmica, 3 (1988), 511-533. 

[33] Rochow H., Analyse und algorithmische Optimierung des MDFS Maximum- 
Matching Algorithmus, diploma thesis, Universitat Bonn (1991). 


52 



[34] Tarjan J. E., Data Structures and Network Algorithms, SIAM (1983). 

[35] Vazirani V. V., A theory of alternating paths and blossoms for proving cor¬ 
rectness of the 0{\/VE) general graph maximum matching algorithm, COM¬ 
BINATORIC A 14 (1994), 71-109. 

[36] Vazirani V. V., A proof of the MV matching algorithm (2014). 

[37] Witzgall C., and Zahn C. T. Jr., Modihcation of Edmonds maximum matching 
algorithm, J. Res. Nat. Bur. Standards, 69 B (1965), 91-98. 


53 



